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Abstract 



At so high temperature (T) that the couphng constant (g) is smah and the masses 
of the particles are neghgible, different scheme has to be apphed in each energy scale in 
the analysis of the quark-gluon plasma (QGP). In the soft energy region (p ~ gT), the 
simple perturbative expansion called the hard thermal loop (HTL) approximation can 
be applied, and that approximation expects the existence of the bosonic and fermionic 
collective excitations called plasmon and plasmino. On the other hand, in the ultrasoft 
energy region (p ~ g'^T), the HTL approximation is inapplicable due to infrared singularity, 
so the question whether there are any excitation modes in that energy region has not been 
studied well. 

In this thesis, we analyze the quark spectrum whose energy is ultrasoft in QGP, using 
the resummed perturbation theory which enables us to successfully regularize the infrared 
singularity. Since the Yukawa model and QED are simpler than QCD but have some 
similarity to QCD, we also work in these models. As a result, we establish the existence 
of a novel fermionic mode in the ultrasoft energy region, and obtain the expressions of 
the pole position and the strength of that mode. We also show that the Ward-Takahashi 
identity is satisfied in the resummed perturbation theory in QED/QCD. 

Furthermore, we derive the linearized and generalized Boltzmann equation for the 
ultrasoft fermion excitations in the Kadanoff-Baym formalism, and show that the resultant 
equation is equivalent to the self-consistent equation in the resummed perturbation theory. 
We also derive the equation which determines the n-point functions with external lines 
for a pair of fermions and {n — 2) bosons with ultrasoft momenta, by considering the non- 
linear response regime using the gauge symmetry. We also derive the Ward-Takahashi 
identity from the conservation law of the electromagnetic current in the Kadanoff-Baym 
formalism. 
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Chapter 1 

Introduction 



1.1 Quark- gluon plasma 

Hadrons are composed of the fermions called quarks and the gauge bosons called gluons, 
and the dynamics of these particles are described by the quantum chromodynamics (QCD). 
The Lagrangian of QCD is as follows: 

C = -\FrF^, + - m,)V'a. {a = 1, 2...N^ - 1) (1.1) 

a=l 

Here Fr = ^^'A''^ - d^AH - gfatcA'^A"^ is the field strength, the quark field, g the 
coupling constant, = + igA'^°'ta the covariant derivative, fabc the structure constant 
of SU{N), nia the current quark mass, respectively. The subscript a in the quark field 
stands for the flavor. In the real world, = 3, Nf = 6 (up, down, charm, strange, top, 
and bottom), and the values for rria are as follows^ [1]: m„ ~ 2 MeV, ~ 6 MeV, rric — 

1.2 GeV, rris ~ 95 MeV, ~ 4.1 GeV, nit ^ 175 GeV. 

Let us briefly review the present theoretical understandings on the phase structure 
of QCD matter. Since the quark number is conserved in each flavor, there exist the 
chemical potential (fii) which are related to the conserved quark number. It is possible to 
analyze the case that the chemical potentials of each flavor are different [2-5], but in the 
following we will review only the case that they have the same value: /ij = /x. We can make 
rather reliable predictions in the case of extremely high temperature or chemical potential, 
because of the asymptotic freedom. At high chemical potential, smallness of g justifies us 
to expect that the system should be a Fermi liquid with a Fermi sphere of the quarks. 
However, there is an attractive interaction between the two quarks in the color anti-triplet 
channel via gluon exchange. This attractive interaction causes an instability of the Fermi 
surface as is well known, and it results in the permanent formation of the Cooper pair 

^The values written here are those of running masses in MS scheme. For it, d, and s quarks, the scale 
is fi; ~ 2GeV while for c, b, and t quarks, the scale is set to that of the mass of each particle. 
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Figure 1.1: schematic figure of a possible phase diagram of QCD. The vertical axis is temperature 
axis while horizontal one is the chemical potential axis. 

composed of quarks leading to the superconducting phase called color superconducting 
(CSC) phase [6-11]. Also in the case that temperature is extremely high, the quarks and 
the gluons are expected to be deconfined [12]. That phase is called quark-gluon plasma 
(QGP) [13]. 

Since g is not so small when temperature and chemical potential are not extremely 
large^, it is difficult to what phase is realized in that region. Nevertheless, the first 
principle calculation using the Monte Carlo simulation, which is called lattice QCD, is 
possible when the chemical potential is zero. The lattice results tell us that the transition 
from the confined phase to the deconfined one is crossover, not the phase transition in the 
usual sense, and that the transition temperature (Tc) is approximately 200 MeV [14-18]. 
Analyses in the small-// region have been tried by using some methods such as the Taylor 
expansion method [19] and the imaginary chemical potential method [20,21]. 

In the case that the chemical potential is finite, the lattice calculation is unreliable at 
present, so we have to use some effective models such as the Nambu-Jona-Lasino (NJL) 
model [22-24] and its improved version [25-27], the chiral random matrix model [28], and 
the strong coupling expansion [29-31]. Many of the analyses [32-35] suggest the following: 
There is a second-order critical point at finite T and fj,. The first order phase transition 
line exist, whose one end is connected with the critical point and the other is on the /i 
axis. Summarizing the suggestions above and the expectations based on the lattice QCD 
and the asymptotic freedom, a schematic figure of a possible phase diagram of QCD is 
drawn in Fig. 1.1. 

^The values for g determined by experiments in some energy scales are as follows [1]: g{l GeV)~ 2.5, 
p(10 GeV)~ 1.5, and ^(100 GeV)~ 1.2. 
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In this thesis, we focus on the analysis of the QGP phase. QGP is considered to be 
reahzed by heavy ion cohision experiment [36,37], which is performed in the Relativis- 
tic Heavy Ion Colhder in Brookhaven national laboratory and Large Hadron Collider in 
CERN. QGP was also realized in the early universe since it is considered that temperature 
was higher than Tc when the age of universe is ~ 10~^ seconds or younger. For this reason, 
the properties of QGP are relevant to the analysis of the early universe. QGP is the only 
experimentally realizable fermion-boson system in which the masses are negligible com- 
pared with temperature. Because of this property, fermionic collective excitations appear 
in QGP, which is the main topic of this thesis. We note that the Yukawa model and QED 
at high temperature have that property, so we also treat these models, which are simpler 
than QCD, in this thesis. Throughout this thesis, the Yukawa model means the model in 
which a fermion (denoted by ip) is coupled with a scalar field (p. The Lagrangian is 

£[0,^,?] = ^(a'^.^)' + V^i(^ + i#)V. (1.2) 
We do not include the possible self-coupling of the scalar fields for simplicity. 

1.2 Hierarchy of energy scales 

In the following, we consider the case that the temperature (T) is so high that the 
current quark mass and the nonperturbative effect are negligible: T ^ mi, Aqqd- Here 
^QCD is a scale parameter, which is of order 200 MeV. The inequality T ^> Aqcb implies 
that 5 <C 1 so that the perturbative expansion in terms of the coupling constant is applica- 
ble. We consider the case that the chemical potential is zero {fj, = 0). Some of rui are very 
large, so it seems that the temperature at which the massless approximation is valid is 
too high to be realized. However, the heavy quarks are almost decoupled from the system 
if their masses are much larger than T, so that approximation is expected to be valid for 
light quarks in that case. One may also have doubt on the validity of the perturbative 
expansion: the coupling constant is not small enough for justification of the perturbative 
expansion in realistic temperature. Nevertheless, there are some discussions [38-44] which 
suggest that the perturbative expansion is at least partially applicable if the temperature 
is higher than ~ 2.bTc. For example, the entropy calculated from approximately self- 
consistent scheme which is motivated by the existence and the properties of the collective 
excitations predicted from the perturbation theory, agrees very well with that calculated 
from the lattice QCD [41-44]. 

In the temperature region described above, the perturbative method suggests that 
the system has multi-energy scale structure [45], which will be reviewed in the following 
subsections and summarized in Fig. 1.2. We note that the structure appears not only in 
QCD but also in the Yukawa model and QED. 

Among the properties of QGP, we focus on the spectrums of the quarks and the gluons. 
These quantity are quite important to elucidate the picture how the quarks and the gluons. 
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Figure 1.2: The energy hierarchy of high temperature QCD. The characteristic quantities, the 
correspondence between the methods in the field theory and the kinetic equations, and the bosonic 
and the fermionic modes obtained from these theories in each scale are shown. The vertical axis 
denote the energy scale. 



which are the basic building blocks of QCD, behave in the system. Naively, that problem 
seems to be trivial in the weak coupling regime: A naive expectation might be that the 
spectrum of the particle seems to be almost the same as that of free particle since g is 
small. We will see that it is not the case when the energy of the particle is ~ gT or g^T 
in the following subsections. 

1.2.1 Hard scale 

Properties of hard particles 

The energy scale T is called hard scale. One example of appearance of this scale in 
the theory is the average inter-particle distance, which can be shown as follows [46]: For 
simplicity, we show in the case of the boson in the Yukawa model. The number density of 
the boson in the free limit is 



) 

d|k||k|2nB(|k|) 



1 

2^ 







= (1.3) 

where nsik^) = l/{exp{k^ /T) — 1) is the Bose-Einstein distribution function, and (^{s) 
is the Riemann zeta function. ((3) is approximately equal to 1.202. We note that the 

1 /3 

volume occupied by one particle is (n^)"^. Since is of order T, the average inter- 
particle distance is of order T~^. 

The de Broglie wavelength of the hard particle, whose momentum is of order T, is also 
of order T~^. This fact suggests that the properties of the hard particle is not so different 
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from those of the free particle because the overlap among the particles is not large enough 
to compensate the smallness of the coupling constant. We show that this suggestion is 
valid in the following. One modification of the dispersion relation of the hard particle is 
the momentum-independent mass called asymptotic thermal mass [42,43,45,47,48], which 
is of order gT. The expressions of the asymptotic thermal masses for the fermions and 
the bosons, in the Yukawa model and QED/QCD, are as follows: 

m} = (1.4) 

ml = (1.5) 

ml = (1.6) 

= (1.7) 

2rp2 

ml = '-^Cj, (1.8) 

Here m'j (ml), ml (m^), and m^ (m^) are the squares of the masses for the fermion (boson) 
in the Yukawa model [47], electron (photon) in QED [48], quark (gluon) in QCD [42,43, 
45,48], respectively, and we have introduced C/ = {N'^ — 1)/{2N). We use g as a coupling 
constant not only in QCD but also in the Yukawa model and QED instead of the standard 
notation (e in the case of QED), to make it clear that the same order counting appears as 
that in QCD in this thesis. 

To demonstrate that the momentum-independent masses appear as a result of the 
perturbative analysis, here we perform an analysis at the one-loop order in the Yukawa 
model. In the present and the next chapter, we employ the real-time formalism in Keldysh 
basis [49,50]. The retarded self-energy of a massless fermion (the boson) S^{p) {11 ^{p)) 
at = is defined as follows: 

^^(P) = , n , ■ ^ n FrTT' (l-^O) 

D^ip) = —o ^n, (1-11) 

where S^{p) {D^{p)) is the retarded fermion (boson) propagator. The pole positions 
of the fermion and the boson are determined by the conditions that the inverse of the 
retarded propagators are zero: 

(/ + ^e- S\p)f - Ipp (l - ^) ' = 0, (1.12) 
p^ + ip^e- n^{p) = Q. (1.13) 
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Note that since /i = and the fermion is massless, Z'^(p) can be expressed in terms of 
the scalar functions S^{p) and as S^{p) = 'yQU^[p) — p • 'yU^ (p) with p = p/|p|. 

Thus we see that the self-energy gives the modification of the pole position due to the 
interaction. 

First, we derive the expression for the asymptotic thermal mass of the boson. The 
boson retarded self-energy at the one-loop order is given by 



-^bare(p) = 



rTr 



'sl{-k)S^{p + k) + S§{-k)Sl{p + k) 



(1.14) 



whose diagrammatic expression is shown in Fig. 1.3. SQ''^{p + k) are the bare propagators 
of the fermion which are defined as 



S^{k) 
S^{k) 



k^ + ik^e' 
1 
2 



npik^Ui^p^ik). 



Here we have introduced the Fermi-Dirac distribution function npik^) 
and the free spectral function p^{k) which is given by 

p°(fc) = 2TTsgn{k°)6{k'^) 



2tt 
2\k\ 



(5(fcO-|k|)-5(fc° + |k|)). 



(1.15) 
(1.16) 

l/(exp(A:0/r) + l) 
(1.17) 



Since the dominant contribution to n^^j.^{p) comes from the region ~ T, as will be 
confirmed later, Tl^^^^ip) ~ g^T"^-, which is smaller than |pp ~ T^. For this reason, 
which satisfies the pole condition (Eq. (1.13)) is expected to be approximately equal to 
±|p|, so p^ is negligible in the calculation of the asymptotic thermal mass. By neglecting 
p^, we get 



d^k 



p\k) 



2g 

9^ 
6 



-7r2 , 



d|k||k|nF(|k| 



(1.18) 



where we have used Eq. (B.l). From this expression and Eq. (1.13), we confirm Eq. (1.5). 

Next we evaluate the asymptotic thermal mass of the fermion. The diagram for the 
retarded fermion self-energy at the one-loop order is drawn in Fig. 1.4. Its expression is 



d^k 

(2^ 



D^{-k)S^{p + k) + D^{-k)S^{p + A;) 



Rf 



(1.19) 
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Figure 1.3: The boson self-energy at the one-loop level. The solid and dashed line are the propa- 
gator of the fermion and the boson, respectively. 



R S 

where Dq ' {—k) are the bare propagators of the scalar boson defined as 

Di{k)=(^l + nB{k')yp'{k). (1.21) 
Inserting Eqs. (1.15), (1.16), (1.20), and (1.21) into (1.19), we obtain 



^Leip)=9'l^yik)[l^l + nsik')^ 



p"^ + 2p ■ k + i{k^ + p^)e 

+ ( -i + ni.(fe°)' 



2 ' J -p2 2p ■ k + i{k^ - p'^)e 

where we have used the on-shell conditions for the bare particles, k"^ = 0, and {k+p)"^ = 
in £)q (— fe) and Sq{p + k). To obtain the expression of the asymptotic thermal mass, we 
only have to calculate {/^, ^barc(p)}' which is easier than calculation of ^^^^j-^ip) when p is 
hard. Let us see this. By using the fact ^barc(^') ~ ^^^i Eq. (1.12) becomes 

p' + iep'-{^,S^,M}=0, (1.23) 

where the terms which are of order g'^T^ have been neglected. Thus we see that the square 
of the asymptotic thermal mass is equal to {^, ^bare(^')}- neglecting as in the boson 
case, {^,^^^^eiP)} '^eads 

J^,X{k), (1.24) 

where 

X{k) = p\k){nF{k^) + UBik"^)). (1.25) 
Note that X{k) is independent of p. By using the formulae Eqs. (B.l) and (B.2), 
{^^'^bire(^')} becomes 

2 roo 

{A^teiP)} = h / ^Ik||k|(n^(|k|) +nB(|k|)) 

-^0 (1.26) 
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Figure 1.4: Diagrammatic representation of the fermion self-energy at one- loop order. The nota- 
tions arc the same as Fig. 1.3. 



From this expression and Eq. (1.23), we confirm Eq. (1.4). 

The other modification of the dispersion relation of the hard particles is the damping 
rate, which comes from the imaginary part of the self-energies. In the Yukawa model, the 
imaginary parts of the self-energies at the leading order have the following forms: 



Imi7^(p)|p, 



-2p°C/(|p|), 



(1.27) 
(1.28) 



where and which will be found to be the damping rate of the fermion and the boson 
at the leading order later, are the real numbers. The leading contribution to and C/ is 
found to come from the two-loop diagrams containing the collision effect among the hard 
particles, not from the one-loop diagrams [47]. The damping rates of the hard particles are 
of order g'^T lii{l/g) (the fermion and the boson in the Yukawa model [47], photon [51-57]) 
or g'^Tliii{l/g) (electron, quark, gluon) [58-62]. The damping rates of the electron, the 
quark, and the gluon at the leading order are known to be independent of momentum [58] . 
The expressions of the damping rates are obtained for the electron (Q), the quark (Cq), 
and the gluon {(g), and only at the leading- log order^ [58-62]: 



C, = 9'TCf 



Co = g^TN 



(1.29) 
(1.30) 
(1.31) 



The asymptotic thermal mass and the damping rate indeed modify the pole position 
of the hard particle. In the case of the Yukawa model, Eqs. (1.12), (1.13), (1.18), (1.26), 
(1.27), and (1.28) lead to 



0\2 



I |2 I 2 
IpI 



2<.(|p|)p°, 



(1.32) 



where i = f or b, and d is the damping rate of each particle. The terms of the asymptotic 
thermal mass and the damping rate, which give the corrections to the dispersion relation 



^In this thesis, we do not regard ln(l/(/) as a large quantity. Therefore ln{l/g) which appears in the 
order estimate of the damping rate is sometimes omitted. 
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and come from the self-energy, are of order g^T'^ and much smaller than the momentum 
term, which is of order T^. This fact is consistent with the expectation that the medium 
correction to the properties of the hard particle is small. In this sense, we can regard the 
fermion and the boson with the hard energy as independent particle excitations. 

Relevance of hard particles to other quantities 

We make a comment on relevance of the hard particle to the physical quantities which 
is not sensitive to the infrared energy region. Whereas the number density of the hard 
particle is of order due to the large phase space as was shown before, that of the fermion 
(boson) whose energy is of order gT is much smaller and of order g^T^ (g'^T^) [45,46], 
as will be shown in the case of the boson in the Yukawa model in the following: First 
we introduce the cutoff parameter A which satisfies A <^ T. From Eq. (1.3), the number 
density of the free particle whose energy is less than A is 



Here we have used nB{x) — T/x for x <^T. This order estimate implies that the dominant 
contribution comes from |k| > A because ns is of order and much larger than nB\k<A- 
By setting A ~ gT, we see that the contribution from |k| < gT is of order g'^T^. This fact 
implies that most of the particles have the hard energy. For this reason, it seems that the 
hard particles determine physical quantities which are not sensitive to the infrared region 
at the leading order. Actually, the leading contribution to some thermodynamic quantities 
(pressure and energy [46,63]) and dynamical quantities (transport coefficients [64-79], self- 
energies of the fermion [80-84] and the boson [80,81,85]), comes from a part in which the 
energy of the particle in the loop integral is hard. We demonstrate this in the case of 
the hard and on-shell boson self-energy. By using Eq. (1.18), the contribution to n^^^^{p) 
from the energy region fc < yl is of order 



Here we have used nF{x) ~ 1/2 for x <^T. This order estimate means that the contri- 
bution from the hard particle, ^ g'^T^, is much larger than that from the particle with 
smaller energy, ~ g^A. 

Nevertheless, we note that there are some quantities which are sensitive to the energy 
region <C T. Thus in analyzing such quantities, it is necessary to take into account the 
contribution from the soft particles even at the leading order. Examples of such quantities 
will be introduced in the next subsection. 




(1.33) 



tA^. 




(1.34) 
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1.2.2 Soft scale 

Properties of soft particles and hard thermal loop approximation 

The energy scale gT is called soft scale. In this scale, the self-energies and the inverses 
of the bare propagators have the same order of the magnitude in Eqs. (1.12) and (1.13), 
as will be shown later. This fact suggests that the medium effects can not be neglected 
even at the weak coupling when the energy scale is soft {p ~ gT). Thus, bosonic and 
fermionic collective excitations can appear in this energy region. There are transverse and 
longitudinal excitations [85] in the bosonic sector, and the latter is called plasmon and 
does not exist in the vacuum. 

To confirm that fact and to see the expression of the dispersion relations and the 
residues of the collective excitations, we here calculate the photon self-energy at the one- 
loop order in QED, adopting the Coulomb gauge. The retarded photon self-energy at the 
one-loop level is given by 



n^Zip) = w'l ^,Tr[rS^{-k)YS^{p + k) + ^^Sl^{-k)YS^ip + k) 



cf I ^,nHk')p'{k)Tr 



r 



^ ( Ik" It + itY:^ ik"^ - i>Y^ 

2k ■ p + p^ 2k ■ p — p"^ 



(1.35) 



where we have neglected the contribution from T-independent part since it can be elimi- 
nated by the renormalization of photon wave function in the vacuum. We note that not 
only the quark but also the gluon contributes to the gluon self-energy in QCD in contrast 
to the case of QED. Since we are focusing on the soft region, we can utilize the useful 
condition p <^k. It is justified since p ~ gT while k T, which was shown in the previous 
subsection. By neglecting the terms which are of order or much smaller than g'^pk, we 
arrive at 



n^Zip)^^9'l^,Mk')Ak) 



^2 



k ■ p \ k ■ p 



(1.36) 



where llIq^ = g'^T^/9. This approximation is called the hard thermal loop (HTL) approxi- 
mation [80,81,86-88]. Here we decompose the photon propagator into the transverse and 
the longitudinal component in the Coulomb gauge: 

^Up) = - (^-^-ZJf^l^ + |p|2+V^L(p)) ' (l-^^) 

where 11 ^^{p) = P'^^^ {p)n^{p)/2 and n^^{p) = —nQQ{p) are the transverse and the 
longitudinal components of the retarded self-energy. Here we omitted the subscript "bare" 
for simplicity. From Eq. (1.36), the retarded self-energy in each component is calculated 
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as 



07 



Q\2 



iPl 



+ 



2|p| 



iPl 



In 



p + IpI 



07 



2|p| 



In 



p + IpI 



(1.38) 
(1.39) 



where we have used Eq. (B.5) in the computation of the transverse component. Thus, the 
dispersion relations in the transverse and longitudinal sectors read 



0\2 



ip' 



P\2 



iPl 



+ 



2|p| 



{P 



0\2 



iPl 



In 



p + IpI 



iPl 



IPI 



■ 3uj, 



07 



1 



In 



p° + IpI 



2|p| 



Ipl 



(1.40) 
(1.41) 



respectively. Since Eqs. (1.40) and (1.41) are invariant for the transformation — )■ —p^, 
we see that the positive energy solution and the negative energy solution are degenerated. 
We analyze only the positive energy solutions from now on. The dispersion relations 
of the excitations in the transverse (a;r(|p|)) and the longitudinal sector (w^dpD) are 
plotted in the left panel of Fig. 1.5. The longitudinal excitation is called plasmon. By 
using p ~ gT, we see that the terms coming from the inverse of the free propagator and 
the terms coming from the self-energy in Eqs. (1.40) and (1.41) have the same order of 
magnitude in contrast to the case of the hard particle. In this sense, the medium effect 
can not be regarded as a weak perturbation in the soft region. The damping rates of the 
two excitations are zero in the HTL approximation, and their leading contribution comes 
from the two-loop diagrams. They are of order g'^T ln{l/ g) [58], so the damping rate of 
the bosonic excitations are much smaller than the excitation energies, which is of order 
gT. The expressions for the residues of the two excitations are given by 



^t(IpI) 
^l(IpI) 



ujt{\p\){ujU\p\) - IpI 



3a;o>2 (IpI) _ (IpI) _ |p 

a;^(|p|)(^i(|p|)-|p|^) 
|p|2(|p|2+3a;2^ 



2\2 • 



^i(lpl)) 



(1.42) 
(1.43) 



The asymptotic forms of the dispersion relations for momenta which satisfies |p| ^ gT 
are as follows: 



4(|P|) 

^i(lpl) 



2 



^^07 



6, ,2 
+ 5"Pl ' 

3, ,2 
+ 5"P" • 



(1.44) 
(1.45) 



Both expressions coincide at |p| =0 because it is impossible to distinguish between trans- 
verse and longitudinal excitations in that case. Prom this expressions, we see that the 
energies of both branches are loq-j at |p| = 0. Thus, we understand the physical meaning 
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p/gT p/gT 

Figure 1.5: The dispersion relations in the bosonic sector and the fermionic sector. Since the 
structure is the same in QED and QCD, we plot the dispersion relations in QED. The vertical axis 
is the energy p", while the horizontal axis is the momentum |p|. The dotted lines denote the light 
cone. Left panel: The dispersion relation in the bosonic sector. The solid line with larger energy 
represents the transverse mode, while the other the longitudinal mode (plasmon). The residue of 
the longitudinal mode becomes exponentially small for |p| ^ gT, so the plot of that mode does 
not represent physical excitation in that momentum region. Right panel: The dispersion relation 
in the fermionic sector. The solid line in the positive energy region corresponds to the normal 
fermion, while that in the negative energy region corresponds to the antiplasmino. Note that since 
we focus on the fermion sector, the antiplasmino appears instead of the plasmino. The residue of 
the antiplasmino becomes exponentially small for |p| 3> gT, so the plot of the antiplasmino does 
not represent physical excitation in that momentum region. 

of ujQ-y: the plasma frequency. On the other hand, for |p| which satisfies gT <ti |p| <^ T, 
the dispersion relations are approximated as 

4(|P|) =^ IpI' + (1-46) 
^l(IpI) ^ IpI (l + 2e"[2|plV(30+i]^ _ (1_47) 

We see that both dispersion relation approaches the light cone as the momentum becomes 
large. However, we note that the residue of the plasmon is exponentially small for large 
momenta: 

ZM) - |PLe-2|Pl^/(3<)-^ (1.48) 

These observations support the expectation that the medium correction is suppressed by 
powers of g when the momentum is much larger than gT. Equation (1.46) tells us that the 
mass of the transverse photon with large momenta is ^y3/2uJQ^ = gT/ which coincides 
the asymptotic mass given in Eq. (1.7). This coincidence is unexpected since the HTL 
approximation is not valid when p ^ T. 

On the other hand, in the fermionic sector, there are also two independent excita- 
tions [82-84]: one is called normal fermion and the other plasmino. The plasmino does 
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not exist in the vacuum as the plasmon does not. Let us perform the calculation of the 
self-energy using the HTL approximation in QED, to show the existence and the proper- 
ties of the fermionic excitations introduced above. The retarded self-energy is expressed 
as 



^LM=i9' I 7^r[D'^^ui-k)S^iP+k)+D^^,{-k)Sgip + k)\j'', (1.49) 
R S 

where DQ^y{—k) is the bare propagators of the photon defined as 

^Om.(^)- + |k|2' (1-50) 

^OM-(fc) = + iPU^)p\k), (1.51) 

where u^^ = (1,0) and P'^y{k) is the projection operator on the transverse direction, 

P'^uik) = g^,igyj{5ij - kikj), (1.52) 

with /c* = /cYl^^l- The diagram corresponding to Eq. (1.49) is shown in Fig. 1.4. Before 
proceeding the calculation, let us show that the term which contains u^u,y in DQ^^{k) does 
not contribute to the result in the HTL approximation. The contribution from that term 
becomes 

-9' j (04/^-H^°)^P°(A:)/. (1.53) 

Here we have used p <^k and neglected the T = contribution. After k integration, this 
expression vanishes. Thus by using the condition A; ^ p as in the bosonic case, Eq. (1.49) 
becomes 

(1.54) 

From this expression, we see that ^\^^^^{p) ~ gT on account of p ~ gT. This order- 
estimate implies that p ~ E^^^^{p), so the computation of {■^, S^^^^{p)} does not give us 
the dispersion relation of the electron in contrast to the case of p ~ T. Therefore we 
proceed the calculation of ^bai-gC^') instead of evaluation of ^bare(p)}- Equation (1.54) 
becomes 

9^ r. ,^ . . 7° - ^ COS • 7 

p'^ — slpl cos 6 

(1.55) 



.(P) = ^9^ I -^,1" [DlA-k)S^{P + k)+ Dl,{-k)Sl{p + k)] r 



^blc(p)=^|^ d|k||k|(n^(|k|)+nB(|k|)) J^^dcose-L^ 



7° 1 + IpI \ P • 7 



■In 



2|p| v^'^ - ipiy ip 



p^ f p^ + \p 



-1 + -win 



2|p| Vp°-IpI 
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In the last line we have used Eqs. (B.l) and (B.2) and introduced 
expression can be rewritten as 



2 

Oe 



Oe 



2|p| 



In 



p° + IpI 



Ipl 



-1 + 



2|p| 



In 



r + IpI 
p'^ — IpI 



. This 

(1.56) 
(1.57) 



By substituting this expression into Eq. (1.10), we obtain 



7 



7 P 



+ ie 



+ 



7° + 7 



^LreiP) 



(1.58) 



From the spinor structure in the numerators of the two terms in the right-hand side, we 
see that these terms are eigenstates of (chirality)/(helicity). The eigenvalue of the first 
term is +1 while that of the second term is —1. In the vacuum, the fermion number of 
the former is +1 while that of the latter is —1. From Eq. (1.58), the dispersion relations 
of the collective excitations which has the fermion number +1 is 



P 



iPl 



Oe 



2|p| 



2+1 



IpI 



In 



p° + IpI 



IPI 



0. 



(1.59) 



The two dispersion relations are obtained as solutions of this equation, which are denoted 
by uj+{p) and —uj-{p). These dispersion relations are plotted in the right panel of Fig. 1.5. 
We note that the excitations whose dispersion relations are w_(p) and — u;+(p) also appear 
in the fermion number —1 sector. The excitation whose dispersion relation is w+(p) is 
called normal fermion, w_(p) the plasmino, and — tj„(p) the antiplasmino, respectively. 
The residue of each branch is 



^±(|P| 



^^(IpI 



IPI 



2< 



Eq. (1.59) can be solved explicitly for |p| which satisfies |p| <ti gT or gT <^ 
In the case of |p| <^ gT, the dispersion relations become 



o^idPl) 



W0( 



1, , 
± - P 



(1.60) 

IpI <■ T. 

(1.61) 



From this expression, we understand the physical meaning of woe^ the counterpart of 
plasma frequency for fermion. The residue becomes 



z±m 



2 3(^0e 



(1.62) 



This expression tells us that both branches have the same strength when they are at rest 
On the other hand, for gT ^ |p| ^ T, the dispersion relations become 



^+(|p|) - IpI + 

a;_(|p|) ~ |p|(l + 2e-2|PPMi-i^ 



(1.63) 
(1.64) 
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and the residue is 



Z+{\P\) 




( 



2|p 




) 



(1.65) 



ZA\P\) 



W -2|p|V-ge-l 

2 



(1.66) 



From these equations, we see that the residue of the plasmino becomes exponentially 
smah as IpI becomes large. We also see that the dispersion relation of the normal fermion 
approaches the light cone and the residue approaches unity, as |p| becomes large. These 
two facts can be understood from the expectation that the medium effect is suppressed 
for T < p, as explained in the previous subsection. Equation (1.63) can be expressed as 
■^+(1?!) — IpP + 2i^0e if neglect higher order term, which implies that the mass of the 
normal fermion with gT S> |p| is equal to the asymptotic thermal mass me at the leading 
order. It is nontrivial because the HTL approximation is not valid when p ^ T. 

We emphasize that the plasmino is a novel excitation which reflects the fact that QGP 
is a fermion-boson system at ultrarelativistic temperature: If each of the mass of the 
quark and the gluon is not negligible, the fermion self-energy would be suppressed and the 
plasmino would not appear. 

Since 2J^{p) ~ gT as was shown before, we note that the fermion retarded self-energy 
and the momentum of the soft particle in Eq. (1.10) have the same order of magnitude 
(~ gT), as in the case of the boson. The damping rates of these excitations are zero in 
the HTL approximation. By evaluating the two-loop diagrams containing the effect of the 
collision, they are found to be of order g'^Tln{l/g) [58]. Thus the damping rates of the 
fermionic excitations are much smaller than the excitation energy cj±(|p|). 

The normal fermion and the plasmino exists not only in /i = case but also in the 
case that fj, is finite [89-91]. It is not straightforward to give an interpretation of the two 
fermionic excitations at finite T since it does not have their classical counterparts, but 
at zero temperature and finite chemical potential case, it is possible to clarify the state 
which forms the plasmino [92]. From that analysis, it was suggested that the state of 
the plasmino is superposition of an antifermion state and the state composed of a boson 
and a hole. The extension of such analysis to the finite temperature case is an interesting 
task. We also note that there is another attempt of interpretation in terms of the idea of 
resonant scattering [93]. 

In addition to collective excitations, the HTL approximation also leads to the Debye 
screening [85]. The screening mass is 




(1.67) 



which is of order gT. 
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Relevance of soft particles to other quantities 

The knowledge of the spectra of the soft excitations is necessary not only to establish 
the picture of the particles, but also to calculate the quantities which are sensitive to the 
soft energy region. Examples include the gluon damping rate at rest. That quantity was 
first calculated by using the bare perturbation expansion and found to be dependent on 
the gauge-fixing [94-98], though that quantity should be gauge-independent [99-101]. The 
solution was given from the following observation: Since properties such as the dispersion 
relation of the soft particles are difi'erent largely from that in the free limit because of 
the medium efi'ect, we have to use the resummed propagators (Eqs. (1.37) and (1.58)) 
which include the information of the result of the HTL approximation instead of the bare 
propagators, when soft particles appear in the loop integral. Vertex function also needs 
to be resummed in such a way. This resummation is called HTL resummation [102-105], 
and it is essentially important to get a sensible result. By using the HTL resummation, 
the sensible and gauge-independent result was obtained [102-105] for the gluon damping 
rate at rest. That method was also applied to the analysis of the quark damping rate at 
zero spatial momentum [106,107]. 

Vlasov equation 

The HTL approximation is a diagrammatic technique for quantities at equilibrium. 
Nevertheless, due to the linear response theory, the HTL results should be reproduced from 
the analysis of the nonequilibrium time evolution at the leading order in the case where 
the system is close to the equilibrium state. Actually, the HTL approximation corresponds 
to the linearized Vlasov equation, which is a collisionless kinetic equation [108-110], as is 
shown in the case of the photon self-energy in the following: The Vlasov equation which 
describes the time evolution of the electron distribution function under the background 
electromagnetic field reads 

vdxnpiKX) = g{E\X) + e'^^v^B^{X))difiF{KX) (1.68) 

where = (1, k) is the four-velocity of a massless particle, ra_F(k, X) the electron distribu- 
tion function at the nonequilibrium state, Ei{X) = —diAQ{X) + dQAi{X) the background 
electric field, Bi[X) = —eijkdxjAi^{X) the background magnetic field with A^[X) being 
the background gauge field, respectively. The left-hand side describes the time evolution 
in the free limit, and is called the drift term. The right-hand side contains the effect of the 
interaction between the electron and the background electromagnetic field, and is called 
the force term. By linearizing Eq. (1.68), we get 

V dx5nF{\^,X) = gw^{X)-^^nF{\k\), (1.69) 

where (5nir(k, X) = ni7'(k, X) — n^^dkl) is the deviation of the electron distribution func- 
tion from the value at equilibrium. We note that the background magnetic field does 
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not contribute in this order. The induced current j^^^^i^) can be expressed in terms of 
(5?ii?(k, X) as foUows: 

f J'^k 

■' ^3 ^ (1.70) 

= ^^V(l^^^"^(''')^-'°^^^^^ + '^^°)' 

where the degeneracy of the spin of electron and the contribution from the positron have 
been taken into account. By using this expression, we get the following expression of the 
polarization tensor by performing the Fourier transformation [X — t- p): 

T= Jo <i|k| V 2 i_i vv + ie) 

Here we have used the relation in the linear response theory, ji^^{p) = U^^'^ {p)Aiy{p), 
and introduced cos0 = p • v. The infinitesimal number +ie comes from the condition 
that the background field is introduced adiabatically: 5nF(k,X) and E(X) vanish when 
— )• —00. This expression yields Eqs. (1.38) and (1.39) after the k integration, which 
means that the linearized Vlasov equation reproduces the photon self-energy in the HTL 
approximation. 

We can also generalize the Vlasov equation to the case where the background field 
is fermionic [108-110]. In this situation, we are to analyze the time evolution of the 
amplitude of the process in which a hard fermion becomes a hard boson and its inverse 
process, instead of the distribution function whose time evolution was analyzed above (see 
Chapter 3). By performing such analyses, the fermion retarded self-energy Eq. (1.55) in 
the HTL approximation, is reproduced. 



1.2.3 Ultrasoft scale 

Resummed perturbation and Boltzmann equation 

The energy scale g^T is called ultrasoft scale. This scale appears in the damping 
rate, or equivalently, the inverse of the relaxation time of the hard particle [58-62], as 
was shown in Sec. 1.2.1. For this reason we expect that the collision effect becomes 
important in the ultrasoft region {p < g^T). To demonstrate that the collision effect is 
not negligible in the ultrasoft region, we consider the linearized Boltzmann equation, which 
is the kinetic equation containing the collision effect, in the relaxation time approximation. 
The situation is the same as in Eq. (1.68): the particle is electron, and the background 
field is the electromagnetic one. That equation reads 

{v-dx + ^ 5n^(k,X) = <7v • E(X)^n^(|k|). (1.72) 
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The second term in the left-hand side is the cohision term in the relaxation time approx- 
imation. Here r is a typical relaxation time for the hard electron, whose inverse has the 
same order of magnitude as the damping rate of the hard electron (Ce), which is of order 
g^Tln(l/^)^ [58-62]. From this equation, we see that the collision term is not negligible 
in the case of ultrasoft region {dx ^ g'^T) since that term and the drift term have the 
same order of magnitude. 

In fact, the Boltzmann equation should be used instead of the Vlasov equation in 
the calculation of the gluon n-point function in this energy region [111-117]. We will 
show in Chapter 3 that the other interaction effects such as the asymptotic thermal mass 
should be taken into account as well as the collision, in the analysis of the ultrasoft 
fermion propagator. We note that the asymptotic thermal mass and a part of the collision 
effect correspond to the real and the imaginary part of Eq. (1.32), respectively. The 
Boltzmann equation in the analysis of the gluon n-point function can be translated into a 
diagrammatic language, and the resultant diagrammatic method is not a simple one-loop 
approximation (HTL approximation), but the resummed perturbation which resums the 
damping rate of the hard particles and sums up the ladder diagrams [111-117]. 

In addition to the gluon n-point function, there are other quantities whose calculation 
need considering the interaction effect among the hard particles. For example, the trans- 
port coefficient is such quantity. That quantity can be obtained from the Kubo formula, 
by taking the zero energy limit of the Green function as follows [63]: The shear viscosity rj, 
the bulk viscosity C, the electrical conductivity a in QED, and the flavor diffusion constant 
Dai3 in QCD are given by 

7? = i- lim i / d^xe'^-''{Mx),7ri„M]), (1-73) 
20 oj-s-o u J 

C = ^ hm - / d'xe^-^\[P{x),Pm), (1.74) 
2 u J 



1 , 1 



a = ^ hm - / d'xe'^-'mx),m]), (1-75) 
D^p = i hm i j d\e^^^\[jf{x),j]m)^-i. (1.76) 

respectively, where P is the local pressure, Tiifn{x) = Ti^ — Si^P the traceless part of 
the energy- momentum tensor with Ti^ being the energy-momentum tensor, ji{x) the 
electromagnetic current, (a;) the flavor current, respectively, and we have introduced 
= 5^0a)- Because the zero energy limit is taken, the energy scale of that quantity 
is much less than g^T. For this reason, the analysis of that quantity also needs includ- 
ing the interaction effect among the hard particles. Accordingly, the transport coefficient 



*For the photon in QED and the particles in the Yukawa model, the mean free path is of order of 
((7^r)^^, so it seems that we do not need to take into account the interaction among the hard particles in 
the ultrasoft region. However, the difference of the thermal masses of the fermion and the boson, which is 
of order g^T^, plays the similar effect to that of the mean free path, as will be shown in Chapter 2, so we 
can not neglect the interaction effect when we analyze the fermion self-energy with ultrasoft momentum. 
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can be calculated either by the Boltzmann equation [72-74] or by the resummed per- 
turbation^ [64-71]. Using the correspondence between that perturbation theory and the 
Boltzmann equation, the resummation scheme in the resummed perturbation theory is 
interpreted with the language of the kinetic theory. 

Suggestion on existence of ultrasoft fermion mode 

When the energy scale is much below g^T, the hydrodynamics works well, and the 
bosonic hydrodynamic modes such as the phonon and the diffusion mode appear. On 
the other hand, the fermionic sector whose energy is of order g^T has not been well 
investigated, due to the difficulty of taking into account of the interaction among the hard 
particles. For this reason, it has not been studied well whether there are any modes in 
this region. 

Nevertheless, there have been some suggestive works for supporting the existence of 
such an ultrasoft fermion mode at finite temperature. Some of them are based on super- 
symmetry, and others are based on analyses using effective model, as will be shown later. 
First we introduce the works related to super symmetry. Historically, the ultrasoft fermion 
mode at finite T was found in supersymmetric models as Nambu-Goldstone fermions called 
goldstino associated with spontaneous breaking of supersymmetry at T 7^ [118], by using 
Ward-Takahashi identity and a diagrammatic technique [119-124]. The complete calcu- 
lation at the leading order was first performed by Lebedev and Smilga [124]. Because 
the supercurrent and the energy-momentum tensor are in the same supermultiplet, the 
goldstino was also regarded as the supersymmetric analogue of the phonon, phonino. Here 
we note that the analysis in [124] was performed in the temperature region T ^ rn/g, 
where m is the bare mass. It implies that their analysis is only valid for gT <^ m, but 
not for m <^ gT for which our analysis in the present paper is concerned. The analysis 
was extended to QCD at so high temperature that the coupling constant is weak [125], 
in which a supersymmetry is still assigned at the vanishing coupling, and hence, the su- 
persymmetry is, needless to say, explicitly broken by the interaction. Thus, there exists 
no exact fermionic zero mode but only a pseudo-zero mode does. Although these analy- 
ses [119-125] are suggestive, it is still obscure whether a genuine ultrasoft fermion mode 
exists when supersymmetry is absent, in particular, at extremely high temperature. 

Here we note that there have been suggestions of the existence of ultrasoft fermion 
mode at finite T even without supersymmetry. It was shown in one-loop calculations [93] 
that when a fermion is coupled with a massive boson with mass m, the spectral function of 
the fermion gets to have a novel peak in the far-low-energy region in addition to the normal 
fermion and the antiplasmino, when T ~ m, irrespective of the type of boson; it means 
that the spectral function of the fermion has a three-peak structure in this temperature 



^That quantity was also computed using the n-particle irreducible formalism [75-79]. All of these 
methods produce the same result in the leading order of the coupling constant. 
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region. It was suggested that such a three-peak structure may persist even at the high 
temperature Hmit in the sense m/T — )■ 0, for the massive vector boson on the basis of a 
gauge-invariant formalism, again, at the one-loop order [126]. Thus, one may expect that 
the novel excitation may exist in the far-infrared region also for a fermion coupled with 
a massless boson, although the one-loop analysis admittedly may not be applicable at 
the ultrasoft momentum region. There are also the works suggesting the existence of the 
ultrasoft fermion mode using the Schwinger-Dyson equation [127-131]; we note, however, 
that it is difficult to keep gauge symmetry in the Schwinger-Dyson approach at finite 
T. The analysis of the quark spectrum around Tc using the NJL model [22-24] was also 
performed [132], and as a result, the existence of the ultrasoft fermion mode was suggested 
due to the coupling between the quark and the mesons. 

Finally let us give a generic argument supporting the existence of an ultrasoft fermion 
mode at finite T on the basis of the symmetry of the self-energy for a massless fermion. 
It was shown that the chiral, parity, and the charge symmetry make the fermion retarded 
propagator have the following structure [133]: 

sV.P) = -iffe^ + ^i?^V (1.77) 



2 \S+(p»,p) S-{p",p} 

Here we have introduced S±{p^,p) = + ie^ |p| — S^{p) ± Z'^(p), and these functions 
satisfy 

5_(/,p) = -(S+(V*,p))*, (1.78) 
S+(p°,0) = 5_(/,0). (1.79) 



By setting |p| = 0, we get 



Eqs. (1.78) and (1.79) imply 



5+(/,0) = -(S+(V*,0))*. (1.81) 

By using 5'+(p°, 0) = p° - 17° (p°, 0), we get 

Rei:°(-/,0) = -Rei:°(/,0), (1.82) 

where p^ is real. This property implies that — KeS^{p^,0), which is the real part of 
the inverse of the retarded fermion propagator, is zero at p'' = if there is no singularity^. 
For this reason, it is suggested that the fermion retarded propagator always have the pole 
around the origin, provided that the imaginary part of the fermion retarded self-energy 



®In the HTL approximation [80,81], the singularity appears at the origin. That is the reason why that 
approximation can not suggest the fermion mode around the origin. 
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is small enough. This argument suggests that the existence of the ultrasoft pole may be 
a universal phenomenon at high temperature in the theory composed of massless fermion 
coupled with a boson, as long as chiral, parity, and charge symmetry exist. We note 
that the attenuation of the pole at the origin in the case that the fermion has finite bare 
mass [134], is consistent with that argument because the finite fermion mass breaks the 
chiral symmetry. 

What we do in the thesis 

In this thesis, we analyze the properties in the ultrasoft region, focusing on the fermion 
spectrum: we calculate the ultrasoft fermion spectrum by using the resummed perturba- 
tion which enables the calculation of the fermion propagator with ultrasoft momentum, 
and show that a novel fermionic mode exists in that energy region. We also obtain the 
expressions for the dispersion relation, the damping rate, and the strength of that mode. 

As explained before, there is an equivalence between the diagrammatic method and the 
kinetic equation such as the equivalence between HTL approximation and Vlasov equa- 
tion, or the resummed perturbation treating the gluon n-point function and Boltzmann 
equation. Because of this equivalence, we expect that the resummed perturbation used in 
the analysis of the ultrasoft fermion is expected to be equivalent to some kinetic equation 
which contains the interaction effect among the hard particles. We derive the generalized 
Boltzmann equation that is equivalent to the basic equation of the resummed perturbation 
theory treating the ultrasoft fermion, from the Kadanoff-Baym equation, which describes 
the time evolution of the nonequilibrium system based on the field theory. 

1.3 Outline of the thesis 

This thesis is organized as follows: Chapter 2 is devoted to the analysis on the spectrum 
of the ultrasoft fermion using the resummed perturbation theory in the Yukawa model 
and QED/QCD. The analysis in the Yukawa model and QED is based on Ref. [135]. As a 
result of the analysis, we show that a novel fermionic mode which we call ultrasoft fermion 
mode exists in that energy region, and obtain the expressions of the pole position and 
the strength of that mode. We also show that the resultant fermion propagator and the 
vertex function satisfy the Ward-Takahashi identity in QED/QCD. In Chapter 3, we derive 
the linearized and generalized Boltzmann equation for ultrasoft fermion excitations, from 
the Kadanoff-Baym equation in a Yukawa model and QED. We show that this equation 
is equivalent to the self-consistent equation in the resummed perturbation theory used 
in the analysis of the ultrasoft fermion spectrum at the leading order. Furthermore, we 
derive the equation that determines the n-point function with external lines for a pair of 
fermions and (n — 2) bosons with ultrasoft momenta in QED. We also showed that the 
Ward-Takahashi identity is satisfied, and that identity can be derived the conservation law 
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of the current. The analysis in that chapter is based on Ref. [136]. Finally we summarize 
this thesis and give the concluding remarks in Chapter 4. In Appendix A, we show that the 
results obtained in the Coulomb gauge in this thesis, also can be obtained in the temporal 
gauge. We write some formulae used in the text in Appendix B. 
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Chapter 2 

Ultrasoft Fermion Mode 



Though there are some suggestion on existence of the ultrasoft mode as is written in 
the previous chapter, It is not a simple task to establish that fermionic mode exist in the ul- 
trasoft region on the complete leading order calculation because of the infrared divergence 
called pinch singularity [64-79,111-117,137,138] that breaks a naive perturbation theory, 
as will be briefly reviewed in the next section. We remark that the same difficulty arises in 
the calculation of transport coefficients [64-79] and the gluon self-energy [111-117] in the 
ultrasoft energy region. Therefore, in this chapter, we analyze the fermion propagator in 
the ultrasoft energy region in Yukawa model and QED/QCD using a similar diagrammatic 
technique in Refs. [124, 125] to regularize the pinch singularity. We shall show that the 
retarded fermion propagator has a pole at = ±|p|/3 — iC (C ~ g'^T\iig~^ for Yukawa 
model and ~ g^T In for QED/QCD) with the residue Z ^ for ultrasoft momentum 
p taking into account the ladder summation. 

This chapter is organized as follows: In Sec. 2.1, we discuss the ultrasoft fermion mode 
in Yukawa model as a simple example without supersymmetry. In Sec. 2.2, we examine 
the ultrasoft fermion mode in QED. We analytically sum up the ladder diagrams giving 
the vertex correction in the leading order, and find the existence of the ultrasoft fermion 
mode as in the Yukawa model. We shall also show that the resultant ultrasoft fermion 
propagator and the vertex satisfy the Ward-Takahashi identity. We perform the analysis 
in QCD using the method which is similar to that in QED in Sec. 2.3. Section 2.4 is 
devoted to a brief summary of this chapter. 

The analysis in Sees. 2.1 and 2.2 is based on Ref. [135]. 

2.1 Yukawa model 

Let us start with the Yukawa model, which is the simplest model to study the ultrasoft 
fermion mode. Generalization to gauge theory will be discussed in Sec. 2.2 and 2.3. We 
calculate the fermion retarded self-energy and obtain the fermion retarded Green function 
with an ultrasoft momentum "p < g^T . We first see that the naive perturbation theory 
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breaks down in this case. Then, we shah show that a use of a dressed propagator gives a 
sensible result in the perturbation theory and that the resulting fermion propagator has 
a new pole in the ultrasoft region. 

2.1.1 Pinch singularity 

The retarded self-energy in the one-loop level is given by Eq. (1.19), and its diagram- 
matic representation is shown in Fig. 1.4. For small p, the self-energy, which is expressed 
as Eq. (1.22) after some manipulations, is reduced to 

where X{k) is defined in Eq. (1.25). This approximation is equivalent to the HTL ap- 
proximation [80,81]. The HTL approximation is, however, only valid for p ~ gT, and 
not applicable in the ultrasoft momentum region. In fact, the retarded self-energy in 
the one-loop level obtained with use of the bare propagators is found to diverge when 
p — 7- 0, since the integrand contains 1/p ■ k. This singularity is called "pinch singular- 
ity" [64-79,111-117,137,138]. 

The origin of this singularity is traced back to the use of the bare propagators because 
the singularity is caused by the fact that the dispersion relations of the fermion and the 
boson are the same and the damping rates are zero in these propagators. For this reason, 
one may suspect that this singularity can be removed by adopting the dressed propagators 
taking into account the asymptotic masses and damping rates of the quasiparticles, as will 
be shown to be the case shortly. 

2.1.2 Resummed perturbation 

Since the leading contribution comes from the hard {k ~ T) internal and almost on- 
shell (/c^ ~ 0) momentum^, we are led to employ the following dressed propagators for the 
fermion and boson: 

- " e-m) + 2K,k«- 
DHk)^[\^n,m) (2.5) 



^Here we note that the case where the internal momenta are soft [k ~ gT) or smaller is not relevant: In 
fact, the HTL-resummed propagators [102-105] should be used for soft momenta. However, the dispersion 
relations of the fermion and the boson obtained from these propagators are different from each other, so 
the pinch singularity will not appear in this case. 
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where the expressions of and mj, are given in Eqs. (1.4) and (1.5). The damping rates 
of the hard particles, Q and Cftj are of order g'^T In g~^. The logarithmic enhancement for 
the damping rate is caused by the soft-fermion exchange, which is analogous to that of the 
hard photon [51-57]. Note that these resummed propagators are the same as those used 
in [125], except for the smallness of the damping rates: We remark that such a smallness 
is not the case in QED/QCD, where the damping rate is anomalously large and of order 
g'^T In ("anomalous damping") [58-62]. 

Using these dressed propagators, we obtain 

for small p, where (5m^ = — m'j = g'^T'^/2A, ( = Q + (h, X{k) = {g"^ / 5m?)X{k)^ and 
= ijp + zC,p)- We have used the modified on-shell condition of the quasi-particles, 
k"^ — rnj + 2i(fk^ = and k"^ — + 2iQW = 0, to obtain the denominator of the 
integrand in Eq. (2.6). We have also neglected rrib, nif, Cfe, and C/ ™ X{k), since the 
leading contribution comes from hard momenta k ^ T. It is worth emphasizing that 
thanks to Sni^ and C,, U^{p) given in Eq. (2.6) does not diverge in the infrared limit, 
p ^ 0. 

Before evaluating Eq. (2.6), we introduce the the following dimensionless value: 

X{k). (2.7) 



(2vr) 



It can be computed as follows: 



/ 







(5m^ J 


(27r)3 






2'K'^5'm? 


Jo ' 


g'T' 




86m^ 





d|k||k|(n^(|k|)+nB(|k|)) (2-^ 



In the last line, we have used Eqs. (B.l) and (B.2). From Eq. (2.8), we see that A is 
of order unity. This value will characterize the strength of residue of the pole for both 
Yukawa model and QED/QCD. 

We expand the self-energy in terms of p^ instead of p^ itself. This is the key point 
of our expansion, which enables us to analytically find the pole of the ultrasoft fermion 
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mode. Then, the leading contribution is 

J4 



X (s|k|7° - k • 7)((p° + iC)s|k| - k • p) (2.9) 

2 foo 

X J (icos0((/ + iC)7° + cos2^p-7). 
By using Eqs. (B.3) and (B.4), we get 

I^''{p)^-^{{p' + iCh'' + vp-j), (2.10) 

with Z = g'^/(8A^7r^) and v = 1/3. Note that the zeroth-order term is absent, which 
implies that there is no mass term. Actually, it is guaranteed by the symmetries, as 
discussed in Sec. 1.2.3. Thus, we obtain the fermion propagator in the ultrasoft region as 

1 



5«(pO,p) 



.0 ^ .,0 , ^ X l2.11) 



Here we have used ^ p) ^ ^ and decomposed the fermion propagator into the fermion 

number +1 and —1 sectors in the second line. These two sectors are symmetric under the 
transformation p o — p and v f-)- — f , so we analyze only the fermion number +1 sector 
in the following. 

From Eq. (2.11), we find a pole at 

pO = -i,|p| - (2.12) 

The dispersion relation of the real part, Rep*^ = — f |p|, is shown in the left panel of Fig. 2.1 
together with the HTL results [82-84] for comparison, where the coupling constant is 
chosen as (7 = 0.1. The imaginary part of the pole reads 

C^g^T\ug-\ (2.13) 

which is much smaller than those of the normal fermion and the antiplasmino [47]. Since 
the real part and the imaginary part of the pole are finite for |p| 7^ 0, this mode is a 
damped oscillation mode. The residue of the pole is evaluated to be 
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Figure 2.1: Left panel: The dispersion relation in the fermionie seetor. In both of the two figures, 
the coupling constant is set to g = 0.1. The vertical axis is the energy while the horizontal axis is 
the momentum |p|. The solid (blue) lines correspond to the normal fermion and the antiplasmino, 
while the bold solid (red) one to the ultrasoft mode. Note that since we focus on the fermion sector, 
the antiplasmino appears instead of the plasmino. The dotted lines denote the light cone. Since 
our analysis on the ultrasoft mode is valid only for |p| <C g^T, the plot for |p| > g^T may not have 
a physical meaning. The residue of the antiplasmino becomes exponentially small for |p| ^ gT, 
so the plot of the antiplasmino does not represent physical excitation for |p| ^ gT, either. Right 
panel: The spectral function in the fermion sector, Eq. (2.15), as a function of energy p*^ at zero 
momentum. 



which means that the mode has only a weak strength in comparison with those of the 
normal fermion and the antiplasmino, whose residues are order of unity. It is worth men- 
tioning that such smallness of the residue is actually compatible with the results in the 
HTL approximation: The sum of the residues of the normal fermion and the antiplasmino 
modes obtained in the HTL approximation is unity and thus the sum rule of the spectral 
function of the fermion is satisfied in the leading order. Therefore, one could have antici- 
pated that the residue of the ultrasoft mode can not be the order of unity but should be 
of higher order. 

The pole given by Eq. (2.12) gives rise to a new peak in the spectral function of the 
fermion as 



TT + V\p\ + iQ 



(2.15) 



which is depicted in the right panel of Fig. 2.1, where |p| is set to zero. Since the ex- 
pression of C, for the Yukawa model is not available in the literature, we simply adopt 
C, = g'^T In g^^ / (2n) in the figure. 
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2.1.3 Suppression of ladder diagrams 

So far, we have considered the one-loop diagram. We need to check that the higher- 
order loops are suppressed by the coupling constant. This task would not be straightfor- 
ward because, Sm'^ ~ g'^T'^ appears in the denominator, as seen in Eq. (2.6), which could 
make invalid the naive loop expansion. The possible diagrams contributing in the leading 
order are ladder diagrams shown in Fig. 2.2 because the pair of the fermion and the boson 
propagators gives a contribution of order 1/(7^, and the vertex gives g^. However, there is 
a special suppression mechanism in the present case with the scalar coupling. 

For example, let us evaluate the first diagram in Fig. 2.2, at small p. The self-energy 
is evaluated to be 



Since there are four vertices and two pairs of the propagators whose momenta are almost 
the same, the formula would apparently yield the factor, X(k)X{l) r-^ x ((5m^)~^ ~ g^. 
One can easily verify that this order estimate would remain the same in any higher-loop 
diagram, so any ladder diagram seems to contribute in the leading order as explained. 
However, this is not the case for Yukawa model with the scalar coupling. An explicit 
evaluation of the numerators of the fermion propagators gives ^(^ — f)f = l/k"^ — ^P, 
which turns out to be of order g^. This is because the internal line is almost on-shell, i.e., 
~ g T^, which comes from the asymptotic masses squared. An analysis shows that 
the same suppression occurs in the higher-order diagrams such as the second diagram in 
Fig. 2.2. Thus, the ladder diagrams giving a vertex correction do not contribute in the 
leading order in the scalar coupling, and hence, the one- loop diagram in Fig. 1.4 with the 
dressed propagators solely suffices to give the self-energy in the leading-order. 

We remark that a similar suppression occurs in the effective three-point-vertex at 
p ~ gT [47] . We also note that this suppression mechanism is quite similar to that found in 
a supersymmetric model for an intermediate temperature region in the sense that gT <^ m 
[124], whereas we are dealing with extremely high-T case. It should be emphasized that 
this suppression of the vertex correction is not the case in QED/QCD, where all the ladder 
diagrams contribute in the leading order and must be summed over [125], as will be shown 
in the following sections. 



Next we explore whether the ultrasoft fermion mode also exists in QED at high T. 
One might expect that the analysis would be done in much the same way as in the Yukawa 
model. It turns out, however, that the analysis is more complicated and involved. It is 
necessary to sum up the contributions from all the ladder diagrams even apart from the 
complicated helicity structure of the photon. In this section, we successfully perform the 




(2.16) 
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Figure 2.2: Some of the ladder class diagrams. The solid and dashed line are the dressed propaga- 
tors of the fermion and the scalar boson, respectively. 



summation of the ladder diagrams in an analytic way, and obtain the fermion propagator 
that is valid in the ultrasoft region. Then we evaluate the pole in the ultrasoft region 
explicitly and examine the properties of the ultrasoft fermion mode in QED. We also 
discuss whether the resummed vertex satisfies the Ward-Takahashi identity. 



2.2.1 Resummed perturbation 

One-loop calculation 

First, we evaluate the contribution from the one- loop diagram. The dressed propaga- 
tors with hard momenta read 



S^ik) . ^--nHnj^ (,._^.;;,^.(,o)2 ^ (2-19) 
0»„(*) = (5 + nsif)) (,,_,„;),;-;^^(,o), . (2.20) 

The expressions of rrie and are given in Eqs. (1.6) and (1.7). The damping rates of 
electron and photon are estimated as Ce ~ [58-62] and ^ g^T In [51-57]. 

Note that (^e is much larger than that in the Yukawa model, which is called "anomalous 
damping" [58-62]. This large electron damping makes the damping rate of the ultrasoft 
mode much larger than in the Yukawa model. Here we have adopted the Coulomb gauge, in 
which the analysis becomes simple thanks to the transversality of the photon propagator. 
We note that the n^u'^/lkl^ term has been omitted in Eqs. (2.18) and (2.20) because that 
term vanishes after the integral, as in Eq. (1.53). 

By using these resummed propagators, the one- loop contribution in the ultrasoft region 
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is evaluated as 

74 



^^ne-iooM = I (^7^ [D%{-k)S''ip + k) + D%{-k)S'{p + k)] Y 
- j {2^Y^ 'l + 2p-k/5m?' 



(2.21) 



where (5?n^ = m^ — mg = — g^T"^ / 12, C = Ce + C^y ~ Ce- Here we have used the same notation 
for 5m? and C, as those in the Yukawa model, although their parametrical expressions are 
different from each other. The factor two in the last line of Eq. (2.21) comes from two 
degrees of freedom of photon polarization. At the one-loop order, we obtain 

^o^ne-loop(p) = -^A2((/ + ZC)/ + • 7) • (2-22) 

We note that this expression has the same structure as that for the Yukawa model; see 
Eq. (2.10). 



Ladder summation 

As already mentioned, the ladder diagrams contribute to the leading-order in contrast 
to the case of the Yukawa model. In this subsection, we sum up all the ladder diagrams, 
and obtain the analytical expressions of the pole position and the residue of the ultrasoft 
mode. 

For this purpose, we introduce the vertex gr^{p, k) defined through the following self- 
energy: 

p r d^k ~ -f^'kPj,(k)r''(p,k) 

J (27r)4 ^ ' l + 2p-k/5n? ^ ' 

Here the vertex contains the contributions from all the ladder diagrams by imposing the 
following self-consistent equation for the vertex function: 
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Equations (2.23) and (2.24) are represented diagrammatically in Fig. 2.3. We have used 
the same approximation as that used in the derivation of Eq. (2.1) for the propagator 
of fermion and photon. We should remark here that this summation scheme using the 
self-consistent equation was first constructed in Ref. [125]. However, we also note that the 
equation has never been solved either analytically nor numerically. In the following, we 
solve this self-consistent equation analytically for small p, and show that the dispersion 
relation does not change from that in the one-loop order even after incorporating all the 
ladder diagrams. 

At small p, Eq. (2.24) reduces to 

^ (p,fe)-7 +y (2vr)4^^^^ + 2(p.A;0/5m2 k ■ k' ' ^^'^^^ 

where we have dropped the term which is proportional to ^ because it only gives higher 
order contribution after being multiplied by ^ in the numerator of Eq. (2.23), because 

Let us solve the self-consistent equation (2.25). We expand the vertex function as 

r^'{p,k) = r^{k) + 6r^'{p,k), (2.26) 

where F^ik) is of order unity and 6r'^{p, k) is of order p/{g^T). 
We first evaluate rQ{k), which can be decomposed as follows: 

r^(k) = 7M(A;) + k''B{k) + uf'Cik), (2.27) 

where A, B, and C are 4x4 matrices. Then the self-consistent equation for /^(^(fc) becomes 

/r/^A-' k'k'^P'^ (k')'yf 

j^X{k' f A{k'), (2.28) 



where B{k) and C{k) in the right hand side vanish due to transversality of the photon 
^i^uik)k'^ = P^( 



propagator: P'^^{k)k'^ = P'^^{k)u^ = 0. By assuming that A{k) is a constant, the integral 



becomes 



(27r)4^^^^^ k-k' ^ 



6m^ J (27r)3 2|k'P ' ' " ' ' k^s\k'\-k-k' 

X (k-7- (k •&')(&' -7))^ 

g2 f d^k' 1, ^A:V-(k'-7)(k-k')„ , 

T^- / tttt^tttt^fC k +nB( k' )7^ — ^ 'A, ^ k 

(^m2 y (27r)3 |k'P "^^1 ""^l (A;0)2 _ (k . k02 

X (k-7 - (k-k')(k'-7))^. 



(2.29) 
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By imposing A;^ = as before, this expression becomes 



9 



5w? J (27r)3 |k 



9 



|k|2 

d|k'||k'|(ni.(|k'|) + nB(|k'|))7^^^^A 



Ik|2 



(2.30) 



where we have used Eqs. (B.l) and (B.2), and dropped 1/: term in the last hne. Then, from 
Eq. (2.28), we find 



A 



1 



B{k) 



1 2A 



7°, C{k)=Q, 



(2.31) 



1 + A' 

where 1 is the unit 4x4 matrix. 

Next we evaluate dr{p, k). From Eq. (2.25), expanding F^ip, k) in terms of 2p-k' / Sm? , 
we find the self-consistent equation for Sr{p, k) as 



6Ff'{p,k) 



-Mk')Pl{k')- 



k-k' 



-AY^^ + sr''ip,k') 



(2.32) 



(2^)4- 

It is not easy to analytically solve Eq. (2.32) directly. So we instead calculate the following 
function: 



snip) 



d^k 



x{k)p^,{kh'^^5r^ip,k). 



(^27r)4"v"/-M^ 
Then Eq. (2.32) leads to the following closed equation 

d'^k 



(2.33) 



snip) 



4A 



+ 



(2vr)4 
d^k' 



X{k) 



(27r 



x{k)pUk)rh 



(iTT) 



k-k' 



k ■ k' 6m? 

X{k')Pj^{k')Y^'5F\p,k') 



{p) - X6n{p). 



(2.34) 



Here we have used Eqs. (2.30) and (2.33) in the second line. The solution to this equation 
is readily found to be dn{p) = —XA'^Uonc-loopip)^ where Eq. (2.31) is used for A. Then, 
the self-energy is evaluated as 



U^ip) 



^Jik)P^^{k)j'^^ (^-^F.^ik) + 6F^ip, k) 

^ "^one-loop 

(p) + 5n{p) 



(2.35) 



1 



(7°(/ + <) + ^^P-7), 



where the residue is 



9 

167r2A2 



(1 + A)^ 



1447r2 ■ 



(2.36) 



2.2. QED 



33 



The pole of the ultrasoft fermion mode has the velocity v = 1/3, damping rate and the 
residue Z. The dispersion of the mode is the same as that in the one-loop level whereas 
the residue is changed. This is our main result for QED. 

2.2.2 Ward-Takahashi identity 

In this subsection, we examine whether the summation scheme, Eqs. (2.23) and (2.24), 
is consistent with the U{1) gauge symmetry. Concretely, we check that our resummed 
vertex function and self-energy satisfy the Ward-Takahashi (WT) identity in the leading 
order of the coupling constant. 

The WT identity reads 

kT^ip, k)=i> + ]/^- S^ip + k)-^ + E^{p). (2.37) 

Since r^j^ijp^k) contains two separated scales k ^ T and p < g'^T, we need to treat them 
carefully. For the hard part, Z!^{p -|- /c) ~ X!^{k) is of order g^T, which is negligible 
compared to ^. In addition, the momentum dependent part U^{p + k) — U^(k) ~ g'^p is 
also negligible compared to E^{p). Therefore the WT identity reduces to 

k,,n\p,k) = ^ + S'^{p). (2.38) 

On the other hand, multiplying Eq. (2.24) by A;^, we have 



k,r^{p,k) = {^ + ft)-^ + 1 ^,x{k')f 



{p + k + k'Y 



d^k' + ^ + , ,/^2^ (2.39) 



1 



where we have dropped the terms of order g^T, and used Eq. (2.23) in the last line. 
This expression coincides with Eq. (2.38), so our self-consistent equation satisfies the WT 
identity. We note that this proof was made without using the expansion in terms oip/ g^T. 
For this reason, Eq. (2.38) is generally valid for p < g'^T, as well as for p <^ g'^T. 

Next, we check whether the explicit solution Eq. (2.31) of the self-consistent equation 
(2.24) at zeroth order in p/{g'^T) satisfies the WT identity to see the consistency with 
the gauge symmetry of the following two conditions adopted to obtain Eq. (2.31): One 
is that terms proportional to ^ in I^(p, A:) is dropped, since they are negligible in the 
self-energy due to = /c^ ~ g'^T'^. The other is that we imposed the on-shell condition 
A;2 = 0, because the internal photon in the self-energy is almost on-shell. Using the same 
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Figure 2.4: Some of the ladder class diagrams which do not appear in the case of the Yukawa 
model and QED. The solid and curly lines are the dressed propagators of the quark and the gluon, 
respectively. 



conditions, we expect that the vertex function satisfies the WT identity in the leading 
order of the couphng constant. In fact, we have for the zeroth order in p 



k^r^ik) ^ltA + k^B{k) ~ 0. 



(2.40) 



Here we have dropped ^ and k"^ in the last equality. 

Finally, we check that the equation determining the vertex function Eq. (2.32), which 
is first order in p, satisfies the WT identity. By multiplying this equation by k^, we obtain 



k^5r^'{p,k) 



d*k' 
W) 



-Xik')PjJk')Y^' 



(2.41) 



where we have dropped ^ as in the previous equation and used Eq. (2.35). Therefore, our 
analytic solution of the self-consistent equation satisfies the WT identity in the leading 
order of the coupling constant. 

We note that without the summation of the ladder diagrams, the WT identity is 
not satisfied when the external momentum of fermion is ultrasoft. By contrast, the ladder 
summation was unnecessary in the Yukawa model, in which the gauge symmetry is absent. 



2.3 QCD 

In this section, we show the existence of the ultrasoft fermion mode, and obtain the 
expressions of the dispersion relation, the damping rate, and the strength in QCD, in a 
similar way to that in QED. The differences between QED and QCD are as follows: First, 
the damping rate of the hard gluon is much larger than that of hard photon. This is 
because the gluon can collides via the soft gluons [61,62] owing to the self-coupling while 
the photon can not [51-57]. Second, new kind of the ladder diagrams shown in Fig. 2.4 
contribute in QCD in addition to the ladder diagrams drawn in Fig. 2.2, again due to the 
self-coupling of the gluons. Despite these differences, we will see that the properties of the 
ultrasoft fermion mode are the same qualitatively in both gauge theories. 

To sum up the ladder diagrams, we use the following self-consistent equation for the 
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^^^^ 




Figure 2.5: Diagrammatic expressions for Eqs. (2.42) and (2.43). The notation is the same as in 
Fig. 2.4. 



quark-gluon vertex function Fa {p., k), which can be derived as in the QED case: 



(27r)4 ^ 'l + 2p-k'/5m'^ "f" ' ^ 'V (k' - k)^ |k'-k|2^ 

X /afec(/"(fc - 2kr + g^'ik + k'r + g^'^ik' - 2kr). 

(2.42) 

Here we have used the Coulomb gauge as in the analysis in QED, and introduced 5m^ = 
ml-ml = g^T'^[N/2A + l/(8iV) + Nf /12] and ( = Cg + Cg- The expression for rrig and 
nig are given by Eqs. (1.8) and (1.9). &iid Q are of order g'^Tln{l/g), which is much 
larger than (^^ due to the self-interaction of the gluon. We note that the presence of the 
third term in the right-hand side, which is absent in the case of QED, is caused by the 
self-coupling of the gluons. The quark self-energy is written in terms of Fa [p., k) as 

^ = y (2^'^^'^ l + 2p.k/5m^ ■ ^'-''^ 

The diagrams for these two equations are drawn in Fig. 2.5. 
At small p, Eq. (2.42) becomes 



(27r)4 ^ 'l + 2p-k'/6m'^ ' \ -2k ■ k' |k'-k|2_ 

X (-2A:"^5-" - 2k'' g^^'^ + g^^'' {k + k'r)fahcW{p, k'), 

(2.44) 

where we have used k"^, k'"^ ~ g'^T^, and the fact that ft and P^a{k) appear in the left of 
the vertex function in Eq. (2.43), and assumed that the vertex function has the same color 
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structure as the bare one [125]: Faip, k) = r^{p, k)ta, where r^[p, k) does not have color 
structure. The third term in the right-hand side in Eq. (2.44) is reduced to 



+ 2p ■ k' /dm? 



7^>' ""n,. (-2fc^V" - 2^^"" + g^'^ik + k') 



-2k ■ k' 



+ i 



1 



{-2k'''g''° - 2k'' g^^ + g^^^^k + k'f 
~ , PLik') 



k' -k|2 

d^k' 



fabctchrpip, k') 



(27r)4 ' 'l + 2p-k'/5m? 



7^^^'— -ffe'^g'''" + feV™ - k'^gP'') 
k ■ k' 



|k'-k|2 



k' - k|2 



g^^ik + k'f 



fabctchrpip, k') 



-X{k') 



+ 



(27r)4 ' 'l + 2p-k'/5m? 
(k'.7-k.7),. 



|k' - k|2fc • k 



k'^'g''^ + k^g^'^ - k^g 
2, , (fc+^OW 



|k'-k|2 



fabctchrp{p,k'). 



(2.45) 



Here we have used the fact that P'^a{k) appear in the left of the vertex function in Eq. (2.43) 
again. By using the fact that ^ appear in the left of the vertex function in Eq. (2.43) again, 
the terms which are proportional to g^''^ are summed up to yield 



-Mk')- 



(27r)4^''"''l + 2p-A;V5?n2 

ffabctctbrP{p,k') 

k-k'\k' -k\^ 



k • 7(A;oA:() - |kV) + k' • 7(^0^0 - |kh 



ko-fo{kok'o - |k'n + k'^-foikok'^ - 



(2.46) 



(27r)4"'"' 'l + 2p- A;V5m2 

ffabctctbrP{p,k') 

^ fc-A;'|k'-k|2 ■ 

This contribution is found to be negligible if we use fc^, k''^ ~ g'^T'^- Thus, Eq. (2.45) 
becomes 



(2.47) 



By using the commutation relation among ta, this expression partially cancels the second 
term in the right-hand side of Eq. (2.44). Then, Eq. (2.44) takes the following form: 



r^^{p,k) = ^p + Cf I ^x{k') 



d^k' ^,,,.k''-f>' + Yk'^' ^' 



k-k' l + 2p-k'/5m? 



P^Jk')rP{p,k'), (2.48) 
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where we have used = Cj. 

We note that this self-consistent equation is the same as that in QED except for Cj in 
the right-hand side: There are no quantitatively novel effects which result from the fact 
that the gauge group is non-abelian. We also note that the same equation was obtained in 
the temporal gauge [125], which is nontrivial if we consider the following facts: In QED, 
the photon propagator in the resummed perturbation scheme is always on-shell, so only 
the transverse component is considered. In that case, apparently there is no difference 
between the Coulomb gauge and the temporal gauge in this scheme. However, in QCD, 
the off-shell gluon appears due to the self-interaction of the gluon. Then, the components 
other than the transverse one appear in the calculation, and the propagator of those 
components are different in both gauges. 

We can solve the self-consistent equation in the same way as in QED. The solution is 
as follows: We expand the vertex function in terms of p/g'^T as in Eq. (2.26). The zeroth 
order solution is 

Here A, B{k), and C{k) are defined by Eq. (2.27), and A is defined by Eq. (2.7). The 
self-consistent equation which determines the vertex function at the first order 5r^{p,k) 
is written as 

sr^iP, k) = c,j ^.xik' f^l-^^;''' ^'pj^ik') [-2^j^AY + 6mp, k')^ . 

(2.50) 

Then 6n{p) defined by Eq. (2.33) becomes 



Sn{p) = "■^ (7°(P° + <) + • 7). (2-51) 

Owing to Eq. (2.43), the retarded quark self-energy is found to be 

^ J (2vr)4^^^^ l + 2p-k/6m^ 

1 



(2.52) 



^-(7°(p° + <) + i;p-7), 



where the expression of the residue Z will be given shortly. We note that this expression 
is the same as that in QED except for the numerical factor. Thus, the expression for the 
pole position of the ultrasoft mode in QCD is the same as in QED, while the residue of 
that mode is not: The residue is 

We can also show that the analytic solution of the self-consistent equation satisfies 
the WT identity in the leading order as in QED, by checking that the counterparts of 
Eqs. (2.40) and (2.41) are satisfied. 
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CHAPTER 2. ULTRASOFT FERMION MODE 



Table 2.1: The expressions of the dispersion relation, the damping rate, and the residue of the 
ultrasoft mode in the Yukawa model, QED, and QCD. Nf and N are set to three. 





Yukawa model 


QED 


QCD 


dispersion relation 




-|p|/3 




damping rate 


Cf + Cb 


Ce 


Cg + Cg 


residue 


5V(72vr2) 




g249/(487r2) 



2.4 Brief summary 

In this chapter, we developed the resummed perturbation theory which was originally 
constructed in Ref. [125] and enables us to successfully regularize the infrared singularity. 
By using this method, we analyzed the quark spectrum whose energy is ultrasoft in QGP. 
Since the Yukawa model and QED are simpler than QCD but have some similarity to 
QCD, we also worked in these models. In QED/QCD, the summation of the ladder 
diagrams had to be done whereas that procedure was unnecessary in the Yukawa model. 
That summation was necessary also from the point of view of the gauge symmetry. As 
a result, we established the existence of a novel fermionic mode in the ultrasoft energy 
region, and obtained the expressions of the pole position and the strength of that mode. 
The expressions for the dispersion relation, the damping rate, and the residue of the 
ultrasoft mode in the Yukawa model, QED, and QCD are summarized in Table 2.1. We 
also showed that the Ward-Takahashi identity is satisfied in the resummed perturbation 
theory in QED/QCD. 
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Chapter 3 

Resummation as Generalized 
Boltzmann Equation 

As is described in Chapter 1, there is a correspondence between some perturbation 
schemes and the kinetic equations: the HTL approximation is equivalent to the Vlasov 
equation [108-110], and the resummed perturbation theory which enables the analysis of 
the ultrasoft gluon is equivalent to the Boltzmann equation [111-117]. Therefore it is 
natural to expect that the resummed perturbation theory which is used in the analysis of 
the ultrasoft fermion in Chapter 2, is equivalent to some kinetic equation. However, we 
note that the equation we will obtain is not the kinetic equation in the usual sense. As 
will be shown later, due to the fact that the excitation we are considering is fermionic, 
not bosonic, that equation describes the time-evolution of the amplitude of the process 
in which the hard fermion becomes the hard boson and its inverse process, not that of 
the distribution function of any particles [108-110]. We call such equation "off-diagonal" 
kinetic equation. 

In this chapter, we derive a off-diagonal and linearized kinetic equation for fermionic 
excitations with an ultrasoft momentum in the Yukawa model and QED, while the Boltz- 
mann equation is derived in the case of bosonic excitations. Our equation is systematically 
derived from the Kadanoff-Baym equation [139], and is equivalent to the self-consistent 
equation in the resummed perturbation theory [125,135] used in the analysis of the fermion 
propagator at the leading order in Chapter 2. The derivation helps us to establish the 
foundation of the resummed perturbation scheme. The kinetic equation will also give 
us the kinetic interpretation of the resummation scheme. Furthermore, we discuss the 
procedure of analyzing the n-point functions (n > 3) not only two-point functions of the 
fermion in QED. 

This chapter is organized as follows: Section 3.1 is devoted to the derivation of the 
generalized and linearized kinetic equation and the discussion on the kinetic interpretation 
of the self-consistent equation in the resummed perturbation theory in the Yukawa model, 
which is the simplest fermion-boson system. In Sec. 3.2, a similar analyses in QED is 
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done in the Coulomb gauge. We also show that the Ward-Takahashi identity is valid in 
this scheme, and that the ?i-point function whose external momenta are ultrasoft can be 
determined by using the gauge symmetry. We briefly summarize this chapter in Sec. 3.3. 

The analysis in this chapter is based on Ref. [136]. We note that a similar analysis to 
that in this chapter can be performed also in QCD [140]. 

3.1 Yukawa model 

In this section, we derive a novel linearized kinetic equation from the Kadanoff-Baym 
equation in the Yukawa model. We will find the vertex correction is negligible, which makes 
the analysis simpler than that in gauge theories. Next, we show that the kinetic equation is 
equivalent to the resummation scheme in the resummed perturbation theory [125,135], and 
discuss the interpretation of the resummation scheme using the correspondence between 
the field theoretical method and the kinetic theory. 

3.1.1 Derivation of the kinetic equation 

Throughout this chapter, we work in the closed-time-path formalism [45,46]. We 
perform the derivation of the kinetic equation in a similar way used in [45, 108-112] by 
applying the gradient expansion to the Kadanoff-Baym equation [139] and taking into 
account the interaction effect among the hard particles in the leading order. 

Let us consider the following situation to analyze the fermionic ultrasoft excitation: 
Before the initial time to, the system is at equilibrium with a temperature T. Then, a 
(anti-) fermionic external source rji^x) (j](x)) and a scalar external source j{x) are switched 
on. As a result, the system becomes nonequilibrium. We will consider the case that j{x) 
and T]{x) vanish and ri{x) is so weak that the system is very close to the equilibrium, i.e., 
the linear response regime. Concretely, we will retain only the terms in the linear order 
of the fermionic average field iP" in the fermionic induced source, which will be introduced 
later. 

Let us consider the generating functional in the closed time formalism [45], 

Z[j,r],rj] = J Vcl^V^Vi^e'^, (3.1) 

with 

S= f d^x[C[(t>,^,ii]-{j4> + i^il + m], (3.2) 
Jc 

where (/> and ij: are the scalar and the fermion fields. The space-time integral is defined as 

[ d^x= [ / (fx, (3.3) 
Jc Jc J 

where j^dx^ is the complex-time integral along the contour C = C~^UC^UC^ in Fig. 3.1. 
We will take to — )• — oo and — )■ oo to factorize out the contribution from the path C^. 
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Rex" 



Figure 3.1: The contour path C in the complex plane. 



The Lagrangian is given in Eq. (1.2). By performing an infinitesimal variation with respect 
lo (j) OY ip in Eq. (3.1), we obtain the following equations of motion: 



-d'^^ix) - gm{x) = j(2;) + jmd(x), 



(3.4) 
(3.5) 



where <P = {(j)) = {ip)) is the expectation value of the scalar (fermion) field, and 
Px['^] = (j^x + ig^ix). Here the expectation value for an operator O is defined as 



(3.6) 



flindix) = g{(j){x)7p{x))c {jmd{x) = g{ip{x)^{x)) c) IS the fermionic (scalar) induced source, 
and the subscript c denotes "connected," i.e.. 



((/>(x)V'(x))c = {(j}{x)ilj{x)) -^{x)<F{x). 



(3.7) 



By differentiating Eq. (3.4) with respect to j{y) and Eq. (3.5) with respect to rj{y), we 
obtain 



i]p^[<P]K{x,y) - gD{x,y)^{x) = i- 



dlK{y, x) - g{^{x){iP{y)i^{x)), + 5(y, x)^^'{x)) = i 



6r]{y) 



Here we have introduced the following propagators: 



D{x,y) = {Tc^{xmy))c = i 
S{x,y) = {Tcip{x)'ilj{y))c = i 



5${x) 



6^{y) 
6f]{x) 



K{x,y) ^ {TcHx)<P{y))c = = 



(3.8) 
(3.9) 

(3.10) 
(3.11) 
(3.12) 
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where Tc means the path ordering on the complex-time path C: 

D{x,y) = ec{x\y'')D>(x,y) + ec{y'',x^)D<{x,y), 

S{x,y)=9c{x\y'^)S>{x,y)-ec{y\x'')S<{x,y), (3.13) 
K{x,y) = ec{x\y'')K>{x,y) + ec{y\x'')K<{x,y), 

with 

D>{x,y)^{(t){x)(t){y))c, (3.14) 

D<{x,y) = {<l){y)<l){x)),, (3.15) 

S>{x,y)^{^{x)^{y))„ (3.16) 

5<(x,y)^(V^(y)V^(x))„ (3.17) 

K>{x,y)^{^{x)m)c, (3.18) 

K<{x,y) = {<l,{y)ij{x)),, (3.19) 

and 9c{x,y) being the step-function along the path C. In the approximations introduced 
later, we can see that K^{x, y) and K^{x, y) coincide, which can be checked by K^{x, y) = 
i9{x^ ,y^)[K^ {x,y) — K^{x,y)\ ~ 0. For this reason, we simply write these two functions 
as K{x, y) from now on. We call K{x, y) off -diagonal propagator" which mixes the 
fermion and boson, while we call D{x^ y) and S{x^ y) ^^diagonal propagators." As will be 
seen in Sec. 3.1.3, in the calculation of the ultrasoft fermion self-energy, the off-diagonal 
propagator is more relevant than the diagonal ones. 

By setting x^ S and y° G C~ in Eqs. (3.8) and (3.9), we obtain 

iP,[.p]K{x,y)-gD<{x,y)^{x) = i^-^J^^, (3.20) 

- d^Kix, y) + g{¥{ymy)i^ix)), + 5<(x, yMv)) = (3-21) 

Here we have interchanged x and y in the second equation. Let us evaluate the right-hand 
side of Eqs. (3.20) and (3.21) using the chain rule: 

Hnd(a^) _ f ^Vindix) 6'F{z) ^ (5r?ind(a;) 6<P{z)\ 



5j{y) Jc V ^^(z) Sj{y) S'^iz) 6j{y)J 

d^z{Uix, z)K{z, y) + S{x, z)D{z, y)), (3.22) 



djindiy) _ f f ^Hz) ^Jmd(y) , '^'f'(^) ^jmd{y) 



dr]{x) Jc \6r]{x) 6(p{z) 6rj{x) S^{z) 

= [ d^z{n{y,z)K{x,z) + S{x,z)E{z,y)). (3.23) 
Jc 



Here we have dropped {5iju\d/ / ^j) and {5jmd/5'^){5^/5fl) since they contain more 
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than one fermionic average field. We have also used 

= (3.25) 

where U (i7) is the fermion (scalar) self-energy [45, 111, 112]. We also introduced the 
off-diagonal self-energy, 

The self-energies are decomposed for arbitrary and on the time path C: 

n{x, y) = ec{x\ y°)7I>(x, y) + x0)77<(x, y), (3.27) 

IJ{x,y) = ec{x\y'')^>{x,y) - Ody^ ,x^)E<{x,y), (3.28) 
S{x,y)=ec{x'',y'^)S>{x,y) + ec{y\x'')S<{x,y). (3.29) 

We have not taken into account contact terms, which is negligible in the leading order. 

Here let us rewrite Eqs. (3.22) and (3.23) in terms of real time integral instead of that 
on the complex-time-path. First we evaluate the diagonal self-energy term: 



d'^zU{x^ z)K{z, y) 

C 



d^zE> (x, z)K{z, y) - / d^zS< (x, z)K{z, y) 



to JxO 

tO-ip 

d'^zU<{x,z)K{z,y) (3.30) 

d'^z{U>{x,z) + U<{x,z))K{z,y) - / d'^zE<{x, z)Kiz,y) 
to Jto 
00 

d'^zU^{x,z)K{z,y). 

-00 

In the last line we have taken — )• —00 and introduced the retarded fermion self-energy 

i:«(x,2/) = i9{x^,y'>)[U>{x,y) + U<{x,y)]. (3.31) 

We used the fact that the term integrated on becomes negligible in this limit [45]. In 
the same way, we get 

d^zn{y,z)K{x,z) c^-i d^zn^{z,y)K{x,z), (3.32) 

Jc J~oo 

where we have introduced the advanced scalar self-energy 

n^{x,y) = -ie{y\x'^)[n>{x,y) - n<{x,y)]. (3.33) 
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Next, we evaluate the off-diagonal self-energy term. The off-diagonal self-energy term 
of Eq. (3.22) becomes 



/ (fzS{x,z)D{z,y) 
Jc 

/ (fzS>{x,z)D<{z,y)+ / d^z~<{x,z)D<{z,y)+ / d^zS<{x,z)D> {z,y) 
Jto JxO Jv° 







d''z{S>{x,z)-S<{x,z))D<{z,y) 

to 

rt°-ip 



+ r d^zS<{x,z){D<{z,y) - D>{z,y)) - [ ' d'zS<{x, z)D> {z,y) 

/oo 
d^z{S^ix, z)D<{z, y) + S<{x, z)D^{z, y)), 
-oo 



(3.34) 



where the advanced boson propagator D^{z,y) = —i6{y^ — z'^)[D^ {z,y) — D^{z,y)] and 
the retarded off-diagonal self-energy S^{x, z) = i6{x^, z'^)[S'^{x, z) — S^{x, z)] have been 
introduced. We stop here and discuss the structure of the off-diagonal self-energy in the 
leading order. 

As we sill see later, we utilize the off-diagonal self-energy in the leading order in the 
linear response regime, which is given by 

S{x,y)=g^S'^{x,y)K{y,x), (3.35) 

where S^{x, y) is the free fermion propagator at equilibrium. The diagrammatic expression 
of this equation is shown in Fig. 3.2. Thus, the components of S are given by 

S^x, y) = y)K{y, x). (3.36) 

Here we perform the Wigner transformation, 

/(fc, X)^ j d^se^^-'f + |, X - I) , (3.37) 

where s = x — y, X = {x + y)/2, and f{x, y) is an arbitrary function. Then, we get 

S^{k,X) = ±g'l 0^S'^{k + k')K{k',X), (3.38) 

with 

S°>{k) = ftp\k){l - npik^)), (3.39) 

S°<{k) = ^p°{k)nF{k^). (3.40) 

As will be seen later, K{k',X) contains 5{k''^). Thus, since we focus on the on-shell 
case A;^ ~ 0, which will be confirmed later, (A: -|- fc')^ ~ 2A; • A' 7^ 0. For this reason. 
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Figure 3.2: The ofF-diagonal self-energy S{k,X) in the leading order. The propagator that is 
composed of the solid line and the dashed line with the black blob is the off-diagonal propagator. 
The other notations are the same as Fig. 1.4. 



S^^{k + k') ~ 0, which implies S^{k,X) ~ 0, so the only nonzero function of the off- 
diagonal self-energy appearing at Wigner- transformed Eq. (3.34) is S^{k,X) ~ S^{k,X). 
Therefore, we drop the second term in Eq. (3.34) because that term becomes negligible 
after the Wigner transformation, and hence the equation becomes 

/• /"OO 

d'^zS{x,z)D{z,y) -i d^zS^{x,z)D<{z,y). (3.41) 

Jc J-oo 

In the same way, we get 

J (fzS{x,z)S{z,y) j (fzS<{x,z)S^{z,y). (3.42) 

Combining Eqs. (3.30), (3.32), (3.41), and (3.42), we get 

/oo 
dS(i;^(x, z)K{z, v) + ^^(x, z)D< {z, y)), 
-oo 

(3.43) 

POO 

-dlK{x,y)+gS<{x,y)^{y)= / dS(77^(z, y)K(x, z) - 5<(x, z)^^(z, y)). (3.44) 



Here we have dropped ^{y){ip{y)tp{x))c because {ip{y)ip{x))c contains more than one 'P. 
Equations (3.43) and (3.44) are the Kadanoff-Baym equations from which the kinetic 
equation is derived. 

After performing the Wigner transformation, Eqs. (3.43) and (3.44) become 

-4 + ^ + ig^{X)^ K{k, X) + igD< {k, X)^{X) 

= i{-U^{k,X)K{k,X) - S^{k,X)D<{k,X)), (3.45) 
(A:^ - ik ■ dx)K{k, X) + gS<{k, X)^{X) 

= n^{k,X)K{k,X) - S<{k,X)S^{k,X). (3.46) 
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Here we have used the following transformation law under the Wigner transformation, 

f{x)g{x, y) ^ f{X)g{k, X) - '-{dxf) ■ {dug) + (3.47) 
f{y)g{x, y) ^ f{X)g{k, X) + '-{dxf) ■ {dtg) + (3.48) 
d^zg{x, z)h{z, y) ^ g{k, X)h{k, X) + Ug, h}p, b. + (3.49) 
where 

{g, h}p. B. = dk9 ■ dxh - dxg ■ dkh (3.50) 

is the Poisson bracket, and neglected higher-order terms that contain dx since we focus 
on the case that the inhomogeneity of the average field is dx ~ g'^T, while a typical 
magnitude of k is of order T. This expansion is called gradient expansion [45, 108-112]. 
We retained the second terms in the left-hand sides of Eqs. (3.45) and (3.46) because the 
first terms, which seem to be the leading terms in the gradient expansion, will cancel out 
in the next manipulation. 

By multiplying Eq. (3.45) by {-i^+^x /'2+ig'P{X)+iU^{k, X)) and adding Eq. (3.46), 
we get 

{2ik ■ dx - 2gl^<P{X) - E^'ik, X)] + n^{k, X))K{k, X) 
= g{^D<{k,X) + S<{k,X))r{k,X). 

Here we have introduced gr{k, X) = g]F{X) + S^{k, X) and neglected higher order terms 
which are of order g'^T'^K and g^T~^r. In the leading order, the coupling dependence in 
D^{k) and S^{k) is negligible, so that D^{k) and S^{k) are replaced by the propagators 
at equilibrium and free limit (g = 0): 

D'^<{k) = p\k)nB{k'^). (3.52) 

S^^{k) is given by Eq. (3.40). We note that though the massless condition fc^ = appears 
in Eqs. (3.52) and (3.40) in the present approximation, k"^ is expected to be of order g^T'^ 
if one takes into account the interaction at equilibrium. For this reason, we will use the 
order estimate k"^ ~ g'^T'^. We also note that K can not be replaced by that at equilibrium 
since K vanishes at equilibrium. 

We see that k'^ terms in the left-hand side of Eq. (3.51) were canceled out and k ■ dx ~ 
g'^T'^ term remains. Thus, we can neglect the terms which are much smaller than g^T^K 
in the calculation of the leading order. Following this line, the diagonal self-energies 
are replaced by those at equilibrium in the leading order, whose diagrams are shown in 
Figs. 1.3 and 1.4: 

{^,E^^'"^\k)] = m), (3.53) 
n^i^'i^k) = ml (3.54) 
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Figure 3.3: The schematic figure of the ofi-diagonal propagator K{x, y) in the Yukawa model in the 
leading order in the linear response regime. The solid (dashed) line with black blob is the resummcd 
fermion (boson) propagator that contains the information on the fcrmion (boson) self-energy S 
(n). The gray blob represents the fcrmionic average field ^. 



as were calculated in Sec. 1.2.1. Note that the imaginary parts of the self-energies ~ 
(7^T^ ln{l/g) and momentum dependence are negligible since they are higher order in the 
coupling constant. We have used the on-shell condition, fc^ ~ 0, which will be verified 
later. 

The same logic as in the diagonal self-energies case justifies substitution of the off- 
diagonal self-energy in the leading order: 

E\k,X) = -g^ j ^^S^R{k + k')K{k',X). (3.55) 

Here S^^{k) is the free fermion retarded propagator at equilibrium, whose expression 
is given in Eq. (1.15). We note that the self-energies can not be neglected in contrast to 
dx ~ gT case^ [108-110], because {^, S}K, UK, T~^S have the same order of magnitude 
as g'^T'^K. 

Using these expressions, Eq. (3.51) becomes 

{2ik ■ dx - 25^(X) + 5m'^)K{k, X) = g^p^{k){nB{k^) + nF{k^))r{k, X), (3.56) 

where 5m'^ = — m^. We note that K{k,X) becomes finite only when /c^ = because 
of d{k'^) in the right-hand side. We also note that liK{k,X) ~ g'^TK{k,X), which is 
confirmed by multiplying Eq. (3.56) by ft from the left. This property makes the vertex 
correction term, (k, X), negligible, which corresponds to the fact that there is no ver- 
tex correction in the analysis using the resummed perturbation theory [135]; see Sec. 2.1.3. 
We also find that the lt<P{X)K{k, X) term is negligible and thus the effect of the average 
field of the scalar vanishes in the present approximation. Thus we get 

i2ik ■ dx + 6m'^)K{k, X) = g^p°{k){nB{k°) + nFik°))^{X). (3.57) 

The schematic figure of K{k, X) is depicted in Fig. 3.3. The solid (dashed) line with the 
blob stands for the resummed fermion (boson) propagator. 



This is because k ■ dxK ~ gT'^K > E}K, UK, T'^E ~ g'^T'^K. 
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3.1.2 Kinetic interpretation 

By introducing the off- diagonal density matris^^ A±{k,X) defined as 



we arrive at the following generalized and linearized kinetic equation from Eq. (3.57): 



It should be noted that this equation is not a usual kinetic equation because ^(k, X) can 
not be interpreted as a distribution function since it is the propagator between different 
particles in the fermionic background. Nevertheless, we call this equation ^^generalized 
kinetic equation!^ because we can obtain the Boltzmann equation if we analyze the time- 
evolution of the diagonal propagator instead of the off-diagonal one [45,111-117]. In fact, 
Eq. (3.59) has the following points which are similar to the Boltzmann equation: 

• The particle is on-shell (/c^ ~ 0). 

• The equation has the common structure as the Boltzmann equation: both have the 
non-interacting part, the interaction part between the hard particle and the average 
ultrasoft field, and the interaction part among the hard particles, which correspond 
to the drift term, the force term, and the collision term in the Boltzmann equation, 
respectively. 

When dx ~ gT, dm? is negligible, and Eq. (3.59) becomes the counterpart of the Vlasov 
equation [45,108-110]. Let us recapitulate the interpretations of each term in Eq. (3.59) 
except for the 5im? term. The first term in the left-hand side describes the time-evolution 
of yl-|-(k, X) in the free limit [g = 0), so this term corresponds to the drift term in the 
Boltzmann equation. On the other hand, the term in the right-hand side expresses the 
effect from the average fermionic field. Hence this term corresponds to the force term in 
the Boltzmann equation. 

Now let us discuss the origin of the dm? term. The origin of the term is 
-{^,I!^{k,X)] + n^(k,X) in Eq. (3.51). The real parts of those terms are of order 
g^T^ while the imaginary parts are g'^Thig~^, so the contribution in the leading order 
comes from the real parts. The difference of the real parts of the diagonal self-energies 
expresses the difference of the dispersion relations of the scalar and the fermion, so we 
call the 5m? term ^^mass difference term." We note that the "mass" here is not the bare 
one but dynamically generated one thorough the interaction among the hard particles. 
(5m^ term has no counterpart in the usual Boltzmann equation, which describes the time- 
evolution of the diagonal propagators, S{k,X) and D{k,X). Therefore we cannot obtain 
kinetic interpretation of that term in the usual sense. 



K{k, X) = 27r6{k'^){e{k°)A+{k, X) + 9{-k' 



M-(-k,X)) 



(3.58) 




(3.59) 
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Here let us see the reason why the mass difference term does not have its counterpart 
in the diagonal case. To this end, we derive the equation that corresponds to Eq. (3.57) 
in the diagonal case. The equation governing the propagator of the fermion is as follows: 



PMS{x, y) + ig{K{y, x))^^'^'^{x) = ^^^^^ {x - y) + i 



5r]{y) 



(3.60) 



Since the second term in the left-hand side contains two 1^, we neglect that term. Now let 
us calculate the right-hand side. We set € and y^ G C~. Since the vertex correction 
term, which contains more than one ^, is negligible, we obtain 



5r]{y) 



d^zX!{x, z)S{z, y) 



c 



to 



d''zE>{x,z)S<{z,y) + / d''zE<{x,z)S<{z,y) 



d''zE<{x,z)S>{z,y) 



to 



d^z{U>{x,z) + U<{x,z))S<{z,y) 



,0 



+ 1 d^zE<{x,z){S<{z,y) + S>{z,y)). 



By taking the limit — )• — oo, we get 
Hnd(a;) 



5r]{y) 



I d\{S''{x,z)S<{z,y) + S<{x,z)S''{z,y)). 



Here we have introduced the advanced fermion propagator, 

5^(x, y) EE -i0(yO, x^)iS>{x, y) + 5<(x, y)). 



(3.61) 



(3.62) 



(3.63) 



By performing the Wigner transformation, we get 



S<ik,X) 



i{U"{k, X)S< {k, X) + U<{k, X)S''{k, X)). 



(3.64) 



The following equation is derived from the conjugate of Eq. (3.64) by using 7*^5^ {x, y)^^ 
S<{y,x): 



S<{k,X) 



4 + ^-^9 Ux) + . {dxHX)) 



(3.65) 



i{S'^{k,X)U<{k,X) + S<{k,X)U'^{k,X)). 
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Here we have introduced the retarded fermion propagator, S^, and the advanced fermion 
self-energy, S^, which are defined as follows: 

5^(x,2/) EE i6ix',y^){S>{x,y) + S<ix,y)), (3.66) 
U^{x,y) = -ie{y\x^){E>{x,y) + E<{x,y)). (3.67) 

By multiplying Eq. (3.64) (Eq. (3.65)) by 

+ ^ + (^cp(x) - idk ■ - iE^'ik, X) (3.68) 

+ ig {^{X) + idk ■ dx^^ - i^^ik, X)^ (3.69) 
from the left (right), we get 

(_A;2 _ ik . dx + 9{2mX) - i^dk ■ {dxHX)) 
+ <I>{X)^x) + E'^{KX)])S<{k,X) 

= -ltU<{k,X)S^{k,X), (3.70) 
S<{k, X){-k^ +ik-thc+ g{2mX) + icfk ■ {dx^{Xm 
-ihc^{X)) + {l^,E^{k,X)}) 

= -S^{k,X)U<{k,X)^. (3.71) 

By subtracting Eq. (3.70) from Eq. (3.71), we get 

{2ik-dx -{lt,U^{k,X)}) S<{k,X) + S<{k,X){lt,U'^{k,X)} 



+ 2g0{X)[S< (k, X), ^] + igid'^x^iXm, dkuS<{k, X)} 
-i5^(X){7^,a^5<(fc,X)} 
: ]/^E<{k,X)S'^{k,X) - S^{k,X)U<{k,X)lt. 



(3.72) 



Here we linearize this equation. By introducing 5S{k,X) = S{k,X) — S^{k) and 
5U{k,X) = U{k,X) - i7(^i)(A:), we arrive at the following equation: 

(2ik ■ dx - {^,u'^^''^^{k) - i7(^'i)^(A:)}) 6S<{k,X) 

- SU^'ik, X)}S'< (k) + S'<{k){ft, 5U^{k, X)} + ig{d-^<P{Xm, d,,S'<{k)} 
= pU<{k,X)S'^^{k) - S°^{k)6U<{k,X)^ 
+ ltU^''''^<{k)6S^{k,X) - 6S^{k,X)u'^^^^<{k)lt. 

(3.73) 

Prom this equation, we see that the corresponding term to the mass difference in the 
diagonal kinetic equation becomes 

-{^, S^^'i^^ik, X)} + {^, i;(°i)^(A;, X)} = -2ilm {^, S^^^^^^ik, X)} (3.74) 



3.1. YUKAWA MODEL 



51 



instead of -{^, U'-'^i^^ik, X)} + i7('=^)^(A;, X), so that the real part is canceled out. For 
this reason, the mass difference term is absent in the diagonal case. 

We note that the term (Eq. (3.74)), which is pure imaginary as a result of the can- 
cellation of the real part, is of order g'^T lii[l/g) and thus negligible since dx ~ g^T. 
We also see that the terms coming from the self-energy have complicated form in the 
diagonal shown in Eq. (3.73), while they are reduced to the simple form, 

-{^,U^'''i^^{k)} + 77('=i)^(fc), in the off-diagonal case in the linear response region. This 
difference comes from the following two facts: 

1. Since we linearize the equation in terms of the deviation from the equilibrium state, 
the terms containing SU and in Eq. (3.73) do not have its counter parts in the 
off-diagonal case [K{k, X) vanishes at equilibrium.] 

2. The structure of the right-hand sides of Eqs. (3.45) and (3.46) after neglecting the 
vertex correction terms are simpler than those of Eqs. (3.64) and (3.65). It is because 



3.1.3 Correspondence between kinetic theory and resummed perturba- 
tion theory 

Here, let us show the relation between our kinetic equation, Eq. (3.57), and the resum- 
mation scheme in the resummed perturbation theory [135] in Chapter. 2. For this purpose, 
we write the fermionic induced source r/ind in terms of the fermion retarded self-energy 
Z!^{p), by using the relation in the linear response theory [45, 108-110] in momentum 
space as 



The expression of r^ind is obtained by using the relation r/ind(-'^) = 9K{X,X) = 
g J d'^k/ (27r)'^ K{k, X) in the present formalism. Thus, from Eq. (3.57), 77ind is 



K<ix,y)=K>{x,y). 



(3.75) 




(3.76) 



By performing the Fourier transformation defined by 




(3.77) 



we obtain 




(3.78) 



Comparing Eq. (3.75) with Eq. (3.78), we obtain the self-energy. 




(3.79) 
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Figure 3.4: The fermionic induced source 77ind(-''^) in the leading order in the Yukawa modcL This 
diagram is obtained by connecting two ends in Fig. 3.3. By truncating we obtain the diagram 
in Fig. 1.4, which expresses S^{p) in the resummed one-loop analysis [135] written in Chapter 2. 



This equation coincides with Eq. (2.6), which is the expression of the retarded fermion 
self-energy [135] except for the absence of the damping rates of the hard particles in 
the denominator. As mentioned in the previous subsection, the damping rates of order 
g'^T\n{l/g) is neglected when the external momentum is of order g^T] one can include 
them by taken into account the imaginary part of Eq. (3.57) if one is interested in the 
damping rate. 

The diagrammatic representation of the fermion retarded self-energy in the present ap- 
proach is the same as that in the resummed perturbation theory [135], which is explained 
as follows: The off-diagonal density matrix yl-|-(k, X), which follows the generalized ki- 
netic equation, is represented by Fig. 3.3. The fermionic induced source 'q\„^{X), shown 
in Fig. 3.4, is diagrammatically obtained by connecting the ends of fermion and boson 
propagators in the right-hand side of Fig. 3.3. This diagram is the resummed one-loop 
diagram [135] appeared in Fig. 1.4 itself except for the fermion average field 



3.2 QED 

In this section, we deal with QED. First we introduce the background field method, 
which is useful to construct the equations for the average fields and the Kadanoff-Baym 
equation in a gauge-covariant form. Next, we derive the generalized kinetic equation in 
the linear response regime adopting the Coulomb gauge fixing, in which the transversal- 
ity of the free photon propagator simplifies the analysis. After the derivation, we show 
the equivalence between the generalized kinetic equation and the resummed perturbation 
theory [125,135], and discuss the interpretation of the terms in the kinetic equation. We 
also check that the Ward-Takahashi identity for the self-energy, consequence of the U{1) 
gauge symmetry, is satisfied in our formalism. Finally we discuss how to compute the 
higher-point-vertex function whose external momenta are all ultrasoft, and make an order 
estimate of it in the weak coupling regime. 
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3.2.1 Background field gauge method 

In the derivation of the average field equation and the Kadanoff-Baym equation in 
QED, it is convenient to formulate them in a covariant form under gauge transformations. 
For this purpose, we employ the background field gauge method [45,95,96,141-145]. In 
this method, the following generating functional is employed: 

Z[j,r],r];A,^,W] = j [Va][V^][Vijy^ , (3.80) 

with 

S= j (fx [C[A^ + a^,^ + ij,W + ^] + £gf - {fa^ + V^r? + 1]^)] , (3.81) 
Jc 

where we dropped the ghost term, which is not coupled with the other fields. A^^ and 
are vector fields, ^ {^) and ip (ip) are (anti-) spinor fields, and is the external current, 
respectively. The Lagrangian of QED has the form. 



C[a,^P,^l^] = -iF'"^[a]F^,[a] + iV'^HV', (3.82) 



where -F^iy[a] = d^a^, — d^a^ is the field strength and D^[a] = + iga^ is the covariant 
derivative. In the background field method, the fields in the Lagrangian are decomposed 
to the classical field, identified as the average fields later, and fluctuations in Eq. (3.80). 
The external sources are chosen to be coupled to a^, ip, and but not to A^, 'F, and 1^. 
We impose the following conditions: 

(a^) = {^) = {^) = 0, (3.83) 

which implies that and ^ {F) can be interpreted as the average parts of the photon 
and (anti) electron field, respectively, and —InZ coincides with the effective action [45, 
95,96,141-145]. 

In the background gauge field method, the gauge-fixing term is chosen to be a func- 
tional of such as 

roP = -A<^M^. (3.84) 

where G[a] is the gauge-fixing function and A is the gauge- fixing parameter. Although 
is fixed by the gauge fixing term, the generating functional, Eq. (3.80), is invariant under 
the background field gauge transformations defined by 

tp'(x) ^ h{x)^{x), A^{x) A^{x) - -h{x)df,h\x), 

9 

ipix) h{x)i;{x), i^{x) lp{x)h\x), a^{x) a^{x), ("^-^^^ 
ri{x) — )• h{x)r]{x), r]{x) fi{x)h^ (x), j^{x) — )• j'^(x), 
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where h{x) = exp[i6{x)]. 

Since the fluctuations covariantly transform under Eq. (3.85), the propagators also 
covariantly transform as 

Df,u{x, y) = {Tca^{x)a^{y))c D^^{x, y), (3.86) 
S{x,y) ^ h{x)S{x,y)hHy), (3.87) 
K^{x, y) = (TcV(x)a^(y)), ^ h{x)K^{x, y). (3.88) 

Also, the Wigner transformed off-diagonal propagator in the leading order of g is 
covariant, which can be confirmed by performing the gradient expansion [45]: 

K^{k,X) ^ h{X)K^{k,X), (3.89) 

which implies that the Kadanoff-Baym equation covariantly transforms with respect to 
the background gauge transformations as will be seen in the next subsection. 

We note that, apart from the covariance with respect to the background field gauge 
transformation, the possible gauge-fixing dependence, which will be confirmed only in the 
Coulomb gauge and temporal gauge in this thesis, should be analyzed. 

3.2.2 Derivation of the kinetic equation 

We work in the Coulomb gauge-fixing condition because this gauge fixing makes the 
analysis simple owing to the transversality of the free photon propagator. The gauge-fixing 
condition is G[a\ = did'' and A — )• oo, which constrains the off-diagonal propagator as 

dlKi{x,y) = Q. (3.90) 

The equations of motion for the average fields are given by 

ilp:,[A\^{x) = r?(x) + 7?ind(2;), (3.91) 

d'^A^ix) - a^aM^(x) - g^{x)Y'^{x) = fix) + ji'^d(^)- (3-92) 

Here we have used Eq. (3.83), and the induced fermionic source and the induced current 
are defined as 

Vmd{x) = g{^{x)il){x))c = g'y^K^{x,x), (3.93) 
jf^j(x) = g{^{x)j>'i;ix)), = gTt{j^'S<{x,x)), (3.94) 

which transform 

Vindix) /i(x)r/ind(x), Ji'ndix) ^ Ji^^ix) (3.95) 

under the background gauge transformations. Therefore, Eqs. (3.91) and (3.92) transform 
covariantly with respect to the background gauge transformation. 
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The equations for the propagators are given by 

/oo 

d^ziU'^ix, z)K^{z, y) + S^{x, z)D<''^^{z, y)), (3.96) 
-oo 

^Q2g^.u _ d^^dnyK,{x,y)+gS<(x,yh^^(y) 

d^zin^^-'iz, y)K,{x, z) - S<{x, y)), (3.97) 



where 11 ^^{x, y) and S^{x, y) are the photon and the off-diagonal self-energies, respectively. 
Here we set x^ G and y^ G C~ . The Wigner transformed equations read 



{S^{k, X)K''{k, X) + ~^{k, X)D<''^'{k, X)), (3.98) 



i-k"^ + ik-dx)K''{k,X) + 



Ko{k,X) + gS<{k,X)j''W{X) 



= n^'"'{k,X)K^{k,X) - S<{k,X)~^^{k,X), (3.99) 

where Eq. (3.90) and the gradient expansion have been used. 

Let us show that is negligible compared with the spatial components. We get the 
following equation by multiplying Eq. (3.98) by 



-ilt + ^ + ig4{X) + iU''{k,X)j (3.100) 

from the left, subtracting Eq. (3.99) from the quantity obtained above, and setting /_f = 0: 
{2ik ■ dx - 2gk ■ A{X) - {lt,U^{k,X)} + {ky)Ko{k,X) - n"^^" {k, X)K^{k, X) 



= gS<{k)f\k,X) 

{kyKo{k,X) - n^'\k,X)Ki{k,X) = gS<{k)r\k,X). 



(3.101) 



Here we have introduced 

gr^'ik, X) = gj''^{X) + "^''{k, X). (3.102) 

Since k^ ~ T, we see that 

K° ~ g'^K\ (3.103) 

Using Eq. (3.103), Eq. (3.99) becomes 

{-k^ + ik ■ dx)K>'{k, X) + kf'k^Koik, X) + gS<{k, X)-i^^{X) 

= ^'^^"'{k,X)K^{k,X) - S<{k,X)S^''{k,X). (3.104) 
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Multiplying Eq. (3.98) by 



-i^ + ^ + ig4iX) + iU''ik,X)] (3.105) 



and Eq. (3.104) by PTi{k) defined below, and subtracting the latter from the former, we 



fit 

obtain 



i-2ik ■ dx + 2gk ■ A{X) + {lt,U^{k,X)})K\k,X) + P^,{k)n^"''{k,X)K^{k,X) 
= -g{^D<''\k,X) + P'"'\k)S<{k,X))r,{k,X), 

(3.106) 



where we have used 



P^i{k)K^{k,X) = -K,{k,X) = K\k,X). (3.107) 



One can show that the background fields and the coupling dependences in the diagonal 
propagators are weak, so that we can replace the electron and the photon propagator, 
which is given as follows, by that in the free limit at equilibrium: 

Dl1;{k) =p\k)nB{k^)P]:Ak). (3.108) 

The diagonal self-energies at on-shell in the leading order are given by 

{^,i7^(°^)(A;)} = ml - 2iCk^, (3.109) 
Pj,(fc)i7^(^'i)"^(fc) = -mlPl^ik). (3.110) 

In contrast to the case of Yukawa model, the damping rate of the electron Q cannot be 
neglected because Qk^ ~ ml, while the photon damping rate of order g^Tln{l/g) can be 
neglected [51-57]. We note that the longitudinal part of the photon self-energy does not 
contribute because the projection operator Pj*(fe) is multiplied. 

The off-diagonal self-energy in the leading order has the following expression, which is 
similar to that in the Yukawa model: 

E^^{k,X) = -g^ j 0^^-^S^^ik + k')j'^KAk',X). (3.111) 

The off-diagonal self-energy in r'^{k, X) has to be retained in the case of QED in contrast 
to the case of the Yukawa model since there is no special suppression mechanism. 
By substituting these expressions in Eq. (3.106), we obtain 

{-2ik ■ Dx - 2iCek^ - 6m^)K'{k,X) = -g^P'^'''{k)p'^{k){nB{k^) + npik^)) 

d^k' k'^-i^ + ^i^y , \ (3-112) 



Here we have used P'^ci{k)Pj°'{k) = —P^u{k) and introduced Sm^ = m^ — rn-g. In this 
gauge-fixing condition, it is apparent that only the transverse component of the thermal 
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Figure 3.5: The diagrammatic representation of self-consistent equation for K^^{x, y) in tlie leading 
order. For simplicity, is not drawn. 



photon contributes to K'^{k,X) because of the projection operator P'^^^[k) appearing in 
the right-hand side of Eq. (3.112). We note that this equation transforms covariantly 
with respect to the background gauge transformation from Eq. (3.89). The diagrammatic 
representation of this equation is shown in Fig. 3.5. 

From Eq. (3.112), we write K^{k,X) in terms of the off-diagonal self-energy for later 

use: 

3.2.3 Kinetic interpretation 

Next, we derive the linearized kinetic equation. Multiplying Eq. (3.112) by 7^ from the 
left and setting A^{X) = 0, we obtain 

{-2ik ■ dx - 2iCek° - 5m^)^{k,X) 

= -g2^p'{k)i7iB{k') + npik^mX) ^3^^^^^ 

We decompose Aj_{\i,X) into positive and negative energy terms as 

Kf'{k,X) = 27r6{k^)[e{k°)A'^{k,X) + e{-k^)A''_{-k,X)], (3.115) 
so that we arrive at the kinetic equation^ from Eq. (3.114): 

2ivdx±^+ 2iCe 4±(k, X) = 2gi,[nB(\k\) + n^(|k|)]l^-(X) 



Ikl 



.s^[nHiki)+nMiki)]pf (.)|: / ^^ ''^'|I^;;[J'^^"" ^^(k^^)^ 

(3.116) 



where we have introduced v'^ = (l,k'). There are two terms that do not appear in the 
Yukawa model analyzed in Sec. 3.1. One is the last term in the right-hand side. The 



^We note that the equation wliicli is similar to Eq. (3.116) was obtained in Ref. [146]. However, the 
authors of Ref. [146] failed to write down all of the leading contributions: their equation does not have the 
third term in the left-hand side and the second term in the right-hand side of Eq. (3.116). 
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Table 3.1: The correspondence between the resummed perturbation theory and the generaHzed 
and hnearizcd kinetic equation. 



Diagrammatic method 


kinetic equation 


thermal mass difference 


mass difference term 


damping rate 


collision term 


ladder diagrams 


correction to force term 



first term in the right-hand side is interpreted as the counterpart of the force term in the 
diagonal case [45,108-110], so the last term in the right-hand side acts like "i/ie correction 
to the force term" at least in the linear response regime. Note that this term mixes the 
positive and negative energy modes, in contrast to the case of the Yukawa model. 

The other is the third term in the left-hand side. This term has a similar form to 
the collision term in the relaxation time approximation of the diagonal case, i.e., pure 
imaginary constant (2z^e) times A±{\<.,X). For this reason, we call this term ^^collision 
term." We note that this term is negligible in the case of the Yukawa model as shown in 
the previous section. 

In the diagonal case [45,111-117], the collision term contains momentum integral for 
the diagonal density matrix. In contrast, such term does not survive in the linearized 
equation in the off-diagonal case. As a result, the collision term in the off-diagonal kinetic 
equation has a similar form to that in the relaxation time approximation. 

We emphasize that the off-diagonal self-energy is not negligible in the off-diagonal 
kinetic equation, while negligible in the diagonal one because the terms containing that 
quantity are higher order in ^. This fact makes the correction to the force term, which is 
absent in the diagonal case, appears in Eq. (3.116). 

As we discussed in Sec. 3.1.2, both of the usual Boltzmann equation and our generalized 
and linearized kinetic equation are composed of the non-interacting part, the interaction 
part between the hard particle and the average ultrasoft field, and the interaction part 
among the hard particles. Which part is the counterpart of the mass difference term, 
the collision term, and the correction to the force term? Because the mass difference 
and the collision term come from the self-energies at equilibrium, they correspond to the 
interaction part among the hard particle. The correction to the force term is a part of the 
interaction part between the hard particle and the average ultrasoft field. 

3.2.4 Correspondence between kinetic theory and resummed perturba- 
tion theory 

Here, let us show the equivalence between Eq. (3.112) and the self-consistent equation 
in the resummed perturbation theory [125, 135] in Chapter 2. To this end, we rewrite 
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Eq. (3.111) in terms of the off-diagonal self-energy using Eq. (3.113): 



X 




f"(A;',X). 



(3.117) 



Here we set = 0. By performing the Fourier transformation, Eq. (3.77), we get 



where r^{k,p) = r^''{k,p)^[p). We note that r^{k,p) is the vertex function [125,135] 
introduced in Chapter 2 whose momenta are hard and ultrasoft. Equation (3.118) is none 
other than Eq. (2.24), which is the self-consistent equation in the resummed perturba- 
tion [125,135]. 

The retarded fermion self-energy is also written in terms of the vertex function: from 
Eq. (3.113), we arrive at 



This expression equals to Eq. (2.23), which is the expression of the fermion retarded self- 
energy in the resummed perturbation theory^ [125, 135]. Thus we see that Eq. (3.112), 
derived in the non-equilibrium state in a linear response regime from the Kadanoff-Baym 
equation, is equivalent to the self-consistent equation in the resummed perturbation theory, 
which is constructed in the thermal equilibrium state. 

Here we discuss the correspondence between each prescription in the resummed per- 
turbation theory [125, 135] and each term in the kinetic equation. As in the Yukawa 
model, the resummation of the thermal mass difference in the resummed perturbation 
theory corresponds to the mass difference term in the kinetic equation. The damping 
rate corresponds to the collision term since both of them contain the damping rate of the 
hard electron, (^e- The ladder summation in the resummed perturbation theory [125,135] 
is caused by the correction to the force term in the kinetic equation. Thus, the ladder 
summation corresponds to the correction to the force term. These interpretations are 
summarized in the Table. 3.1. 

^Equation (3.118) is analytically solved for the energy region \po + iCe\ "C g^T in the previous chap- 
ter [135]. 




(3.118) 



X 



l^'Pl{k')p^{k')[nB{k'^) + n^(fc^O)]^"(fc^p) 
2k' ■p + 2iCek'° + 6m^ 




(3.119) 
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3.2.5 n-point functions {n > 3) 

The fermionic induced source rj[^d generates the n-point functions with n > 3, not 
only the fermion self-energy. In this subsection, we derive the self-consistent equation 
determining the ra-point function whose external lines consist of two fermions (iF) and 
(n — 2) bosons (A^) with ultrasoft external momenta, and make an order estimate of the 
quantity. For example in the case of n = 3, we obtain the correction to the bare three-point 
function, gji,6^'^\p -q-r), from rn^^ [45,86-88, 108-110]: 



6^{q)5A^{r 



4 , (3-120) 

= 'J WfSmSA.ir)'^^''^^^- 

Here we have expanded K^{k,X) around A^^ = 0: 

K^'{k, X) = Rf'ik, X)a=o + SKf^ik, X) + 0{A^^), (3.121) 

where K^^{k, X)^=o contains one ^ while 5K^{k, X) contains one ^ and one A^. 

6K^{k,X) can be obtained by the following way. Collecting terms that contain one 
A^' in Eq. (3.112), we obtain 

{-2ik ■ dx - 2<efc° - 5m^)5K\k, X) + 2gk ■ A{X)K'' {k, X)a=o 

/Mu t,a i^au (3.122) 
^.^^T^^^^^(k',x). 

Since K^^{k, X)^=o is determined by setting A^^ = in Eq. (3.112), this equation is closed 
and 5K^^{k, X) can be determined. 

Let us estimate the order of 5Kt'{k,X). From Eq. (3.112) and K^^q ~ g-'^T'^'F, we 
find 

dK^ ~ g~^T~^A^'K'^_n 

o . (3-123) 
- g-^T-^FA'^. 

Therefore, the vertex correction is estimated as 

g6r''r^g~\ (3.124) 

which is much larger than the bare vertex, gj^ ~ g. Similar order estimate for the n-point 
function with n > 3 can be done with the same procedure; as a result, we find that the 
order of the n-point-vertex function is . 
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3.2.6 Ward-Takahashi identity 

We show that the off-diagonal self-energy given in Eq. (3.111) satisfies the Ward- 
Takahashi (WT) identity. From Eq. (3.111), we get 

Here we have used l/t'K^{k', X) = 0, which can be confirmed by multiplying Eq. (3.112) by 
It from the left. This equation generates the WT identity (Eq. (2.38)) derived in Chapter 2 
by setting A'^ = and differentiating with respect to 1^. The WT identity implies that the 
vertex correction is not negligible because the identity relates the vertex correction to the 
fermion self-energy, which is much larger than the inverse of the free fermion propagator 
with an ultrasoft momentum. In the Yukawa model, the WT identity associated with 
gauge symmetries is absent from the outset, so the smallness of the vertex correction is 
not in contradiction with any identity. 

Equation (3.125) can be derived from the conservation law of the induced current, 

-igiVind^ - ^nmd){x) = d^j^^^ix). (3.126) 

By differentiating Eq. (3.126) with respect to ^{y), we obtain 

- g (^^^^^^{x) - 5^{x^ - 2/0)5(3) (x - y)w(x)) = d^S^^iy^x). (3.127) 

Here 5'~' {x — y) is the delta function defined along the contour C. By multiplying this 
equation by J (i^sexp(iA; • s), we get 

- g [ d'se^'-^ (^JJ^^^x) - 5{x^ - y')5^^\^ - y)mUx)] 

J V my) J (3.128) 

= j d'se'''-'d;S^^{y,x). 

Here we have set x^, y^ £ C"*" and neglected the sub- leading terms. The first term in 
the left-hand side has the same order of magnitude as the hard fermion self-energy S{k) 
times ^{X), so that term is negligible. Thus the left-hand side becomes giji^diX). The 
right-hand side becomes 

/O /-oo 
d^se^^-'dl^S>^'{y,x)+ / d^se^^-'dl^S<^'{y,x) 
-00 JQ 

= r dSe*'=-^5^S>^(y, x) + r d'se^'^-'d'^iSO^iy, x) - S>'^{y, x)) ^^-^^O) 

J-oo Jo 

= -^k^S>^'i-k, X) - k^S^^{-k, X). 
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We see that the first term in the last line is negligible because of the on-shell condition. 
Thus we obtain 

gVmdiX) = -k^S'^^'i-k, X), (3.130) 
if we remember that S^{k,X) ~ S^{k,X). This equation is nothing but Eq. (3.125). 

3.3 Brief summary 

In this chapter, we derived the linearized off-diagonal Boltzmann equation for the 
ultrasoft fermion average field in the Kadanoff-Baym formalism, and showed that the re- 
sultant equation is equivalent to the self-consistent equation in the resummed perturbation 
theory at the leading order of the perturbative expansion, while it has been known that 
the Vlasov equation is equivalent to the basic equation of the HTL approximation. This 
equivalence shows that the resummation scheme of that perturbation theory contains the 
interaction effect among the particles beyond the mean field approximation. We note 
that the off-diagonal Boltzmann equation contains not only the collision term but also 
the thermal mass terms and the off-diagonal self-energy term, in contrast to the usual 
Boltzmann equation. We also derived the equation which determines the n-point func- 
tions with external lines for a pair of fermions and {n — 2) bosons with ultrasoft momenta, 
by considering the non-linear response regime using the gauge symmetry. We also derived 
the Ward-Takahashi identity from the conservation law of the electromagnetic current in 
the Kadanoff-Baym formalism. 
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Summary and Outlook 



In this thesis, we developed the resummed perturbation theory which takes into ac- 
count the separation of the scale and enables us to successfully regularize the infrared 
singularity. By using that method, we analyzed the spectral properties of the ultrasoft 
fermion, including its existence itself, in the Yukawa model and QED/QCD. The proce- 
dure of that method consists of the resummation of the asymptotic thermal masses. In 
QED/QCD, the summation of the ladder diagrams is also necessary. As a result of the 
analysis, we established the existence of the novel fermionic mode in that energy region, 
and obtain the expressions of the pole position and the residue of that mode, which are 
summarized in Table. 2.1. We also showed that the resultant fermion propagator and 
the vertex function satisfy the Ward-Takahashi identity in QED/QCD. We also derived 
the linearized and generalized Boltzmann equation for ultrasoft fermion excitations, which 
describe the time-evolution of the amplitude of the process of changing the fermion to the 
boson, from the Kadanoff-Baym equation in a Yukawa model and QED. We showed that 
this equation is equivalent to the self-consistent equation in the resummed perturbation 
theory used in the analysis of the ultrasoft fermion spectrum at the leading order. It 
helps us to establish the foundation of the resummed perturbation theory. We derived 
the Ward-Takahashi identity from the conservation law of the electromagnetic current in 
QED. Furthermore, we derived the equation that determines the n-point function with 
external lines for a pair of fermions and (n — 2) bosons with ultrasoft momenta in QED, 
by using U{1) gauge symmetry. 

We analyzed the finite temperature and zero chemical potential (fx) case in this thesis, 
so it is quite natural to ask how our result on the ultrasoft fermion mode is modified in 
other region. Such analysis is already done by the author and his collaborator [147], and 
we hope to report its results elsewhere. We just describe the analysis we performed and its 
results in the following: We constructed the resummed perturbation theory in high-// and 
T = case, and analyzed the fermion spectrum with the ultrasoft momentum (< g'^fi) in 
the Yukawa model. As a result, we found that the ultrasoft fermion mode does not exist in 
this case. This result is consistent with the absence of the charge symmetry in high-// case, 
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because the charge symmetry is essentially important for generating the ultrasoft fermion 
mode, as discussed in Sec. 1.2.3. We also investigate how large /i kills the ultrasoft fermion 
mode, and found that the maximum chemical potential in which the mode persists is of 
order T. 

The resummed perturbation scheme used in this thesis can be applied to other quanti- 
ties in the low energy region, in general multi-component system. We briefly discuss some 
examples in the following: 

• As we wrote in Chapter 1, the computation of the transport coefficient also needs 
taking into account the interaction effect among the hard particles. Usually, only the 
collision effect need to be resummed. However, our analysis in this thesis suggests 
that the difference of the dispersion relations of the particles also needs to be taken 
into account, when two kinds of particle appear in the loop integral. It is the case 
in the computation of the flavor diffusion constant [72]: the quarks with different 
flavor appear in the integral. Since the current quark masses are different in each 
flavor, it is quite interesting to investigate the effect of the difference of the current 
quark mass by performing the analysis which takes into account the current quark 
mass. 

• The perturbative analysis of the neutrino spectrum [148-150] around the critical 
temperature of the electroweak transition, which is expected to be realized at the 
early universe, needs the resummation used in this thesis since the bare masses of the 
weak bosons becomes negligible compared with the temperature. It is an interesting 
task to perform such analysis, and investigate the effect of the fermion spectrum in 
the ultrasoft region, on the properties of the early universe. 

• Another example is the supersymmetric system in both of the relativistic and non- 
relativistic system. In such system, the SUSY is spontaneously broken by the fi- 
nite temperature/density effect [118,151], and as a result, the fermionic zero mode, 
goldstino, appears [119-124,151]. As is described in Sec. 1.2.3, the analysis of the 
goldstino, which uses the resummation scheme, was performed in relativistic sys- 
tem [119-124]. 

Though the SUSY has not been observed as a fundamental symmetry in the real 
world, it was suggested to simulate the supersymmetric system using the cold atom 
system [151-158], which is nonrelativistic. Such system also needs the resummation 
scheme [158], and it may be surprising since the power counting rule in the non- 
relativistic system is different from that in the relativistic system. We note that, 
however, there is an important difference between the relativistic system and the 
nonrelativistic one: Since the energy-momentum tensor and the supercurrent are 
not in the same supermultiplet in the nonrelativistic system, the goldstino can not 
be interpreted as the phonino. The dynamical analysis on the goldstino in the non- 
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relativistic system was performed [158], but without the Bose-Einstein condensation 
(BEC). The analysis in the system with BEC is quite attractive task. 
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Appendix A 



Kadanoff-Baym Equation in 
Temporal Gauge 



In this thesis, the calculation was performed in the Coulomb gauge. However, since 
the resummed perturbation theory was first proposed in the temporal gauge [125], we need 
to check that the results in this thesis also can be obtained in the temporal gauge. To 
this end, we show that the equation determining K^^ in the temporal gauge is the same as 
that in the Coulomb gauge, Eq. (3.112), in QED in this Appendix. Showing that equation 
can be obtained also in the temporal gauge suffices to show that the expression of the 
properties of the ultrasoft fermion mode obtained in the Coulomb gauge do not change in 
the temporal gauge, because that equation is equivalent to the self-consistent equation in 
the resummed perturbation theory, which was used to analyze the ultrasoft fermion mode. 

The gauge- fixing condition in the temporal gauge is G[a] = and A — )• oo. This 
condition is equivalent to the constraint = 0. Because of this constraint, we have 



K\x,y) = D'>''ix,y)=0. 



(A.l) 



The equations governing are 



= -i{E^{k,X)K\k,X) + Sf{k,X)D<^\k,X)), 



(A.2) 



{-k^ + ik- dx)K'\k,X) + (k 



idx 



2 



) 



( 




(A.3) 



+ gS<{k,X)^'^{X) 

n^'\k,X)Kj{k,X) - S<{k,X)S^'{k,X). 
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From these equations, we obtain 

{-2ik ■ dx + 2gk • A{X) + {ft, S^{k, X)})K'{k, X) 

-{k-'^y(k- '^yK,ik,X) + n^^^ik,X)K,ik,X) (A.4) 

= -{ltD<^\k,X) +6'^S<{k,X))rj{k,X). 

Here let us evaluate k'^Ki{k, X), which is the longitudinal component of K^(k,X). By 
multiplying Eq. (A.4) by ki, we get 

|k|2 (^k - '-^y Kj{k, X) + hn^'^k, X)Kj{k, X) = -S<{k, X)hn{k, X). (A.5) 

Here we have neglected the terms that are of order g'^T^k^Ki(k,X). We see that 
k^Ki{k,X) ~ g^TKi{k, X) and thus the longitudinal component of Ki{k,X), k^Ki{k, X), 
is negligible compared with the transverse component of Ki(k,X). We note that = 0, 
which is the result of the gauge-fixing condition, and k'^Ki = are valid also in the Coulomb 
gauge in the leading order. Furthermore, also the free photon propagator at equilibrium 
is the same as that in the Coulomb gauge. Thus, we can obtain Eq. (3.112) in the same 
way as in the Sec. 3.2.2. 
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Appendix B 

Formulae 



In this appendix, we write some formulae used in the text. 



/■OO ZrpZ 

d|k||k|n^(|k|) = ^, (B.l) 

d|k||k|nB(|k|) = -^, (B.2) 

d|k||k|V(|k|) = ^^, (B.3) 

d|k||k|3nB(|k|) = ^. (B.4) 



^2 qS^ 

dx = -^^ — —^ + 2a{x-a) + a'^\n{x-a). (B.5) 



X — a 
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